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Abstract 


Low  cycle  fatigue  (LCF)  of  laminate  composite  structures  used  in  Army  applications  is 
assessed  to  identify  the  key  physical  phenomena  occurring  during  LCF  processes  and  to 
determine  their  main  characteristics.  Special  attention  is  given  to  the  LCF  conditions  inherent  in 
Army  structures  (i.e.,  high  cyclic  or  pulse  loads  reaching  up  to  90%  of  the  ultimate  strength  in  a 
fraction  of  a  second).  A  summary  of  fatigue-related  issues  in  laminate  composites  employed  in 
Army  land  combat  systems  is  presented.  Analysis  indicates  that  finite  strain  rate  effects  are 
important  under  LCF  conditions  and  the  pulse  vibration  fatigue  (PVF).  Fatigue  damage 
mechanisms,  evolution  patterns  of  damage,  and  damage  accumulation  processes  are  singled  out 
and  thoroughly  analyzed  as  the  key  mechanical  phenomena  contributing  to  the  changes  in  the 
material  damage  state  and  the  property  degradation  under  fatigue  conditions.  Possible 
correlation  between  ballistic  and  LCF  performance  is  discussed.  Various  models  for  damage 
accumulation  and  fatigue  life  predictions  are  reviewed.  Recommendations  for  fundamental 
research  in  the  areas  relevant  to  the  LCF  of  composite  structures  are  included  to  establish  a 
conceptual  framework  for  the  U.S.  Army  Research  Laboratory  (ARL)  LCF  Program. 
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1.  Introduction 


Advances  in  composite  manufacturing  technology  lead  to  significant  improvements  in  the 
mechanical  properties  of  composite  materials  used  in  defense  industries.  Better  control  of 
processing  defects,  reduction  of  void  content,  and  optimization  of  interphase  properties  improve 
structural  performance  of  composites  and  yield  higher  strength-to-weight  ratios,  stiffness,  and 
fracture  toughness.  However,  multifunctional  hybrid  composite  materials  with  further  enhanced 
properties  must  be  developed  to  meet  the  unprecedented  needs  of  Army  land  combat  systems 
for  speed  and  mobility  without  sacrificing  ballistic  protection  (Bums  1998). 

The  newly  developing  hybrid  composite  structures  consist  of  a  number  of  dissimilar 
materials,  which  themselves  represent  complex  composite  systems.  The  mechanical  properties 
of  novel  hybrid  composites  are  optimized  by  designing  and  tailoring  various  types  of  interphases 
at  the  molecular  and  nanostructural  level.  The  properties  of  polymers  adjacent  to  solid  surfaces 
can  be  significantly  perturbed,  which  can  result  in  dramatic  changes  in  bulk  composite  material 
behavior  (Palmese  and  McCullough  1994;  Fink  and  McCullough  1999;  VanLandingham  et  al. 
1999).  New  materials  with  functionally  graded  properties  can  be  used  to  achieve  optimal 
structural  performance  of  the  designed  hybrid  systems  (Suresh  and  Mortensen  1998).  It  is 
imperative  that  such  advanced  composite  structures  possess  superior  mechanical  properties  over 
a  long  period  of  service  life  in  harsh  environmental  conditions.  The  full  potential  of  composites 
as  advanced  materials  can  be  realized  only  if  the  deterioration  of  long-term  material  properties  can 
be  properly  understood  and  controlled.  Poor  fatigue  performance  may  significantly  reduce  the 
weight  advantage  of  many  composites  (Reifsnider  1991;  Talreja  1993).  This  is  especially 
important  for  the  Army  systems  subjected  to  the  low  cycle  fatigue  (LCF)  conditions  (i.e.,  high 
cyclic  or  pulse  loads  reaching  up  to  90-95%  of  the  ultimate  strength  in  a  fraction  of  a  second). 

In  addition  to  novel  design  strategies,  the  modeling  of  composite  material  structures  must  be 
enhanced  to  enable  cost-effective  design  practices.  New  modeling  capabilities  rooted  into 
fundamental  research  (Talreja  1987,  1996;  Reifsnider  1991;  Case  and  Reifsnider  1998)  are 
especially  needed  for  simulation  of  the  high-load  LCF  environment  inherent  in  many  Army  land 
combat  systems.  To  extend  the  life  cycle  of  new  hybrid  structures,  the  deterioration  of 
composite  properties  such  as  strength,  toughness,  durability,  and  impact  resistance  must  be 
understood  and  controlled.  The  advanced  modeling  capabilities,  when  integrated  with 
appropriate  design  philosophies  and  methodologies,  will  result  in  significant  performance 
improvements  of  vehicles,  ordnance,  and  other  Army  polymer  composite  structures. 
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The  purpose  of  this  report  is  to  assess  the  material  issues  affecting  the  fatigue  behavior  of 
various  composite  structures  used  in  Army  applications  and  to  identify  the  key  physical 
phenomena  present  in  fatigue  processes  and  their  main  characteristics.  Special  attention  is  given 
to  the  effects  of  LCF  conditions  inherent  in  Army  systems.  A  summary  of  fatigue-related  issues 
in  laminate  composites  used  in  Army  applications  is  presented  in  section  2.  Following  Talreja 
(1987)  and  Reifsnider  (1991),  the  fatigue  damage  mechanisms  and  the  kinetics  of  fatigue  damage 
accumulation  processes  are  singled  out  as  the  key  mechanical  processes  controlling  the  property 
degradation  rates  under  fatigue  conditions  (section  3).  Various  models  for  damage  accumulation 
and  fatigue  life  predictions  are  reviewed  in  sections  4  and  5,  respectively.  General 
recommendations  for  fundamental  research  in  the  areas  relevant  to  the  LCF  of  composite 
structures  are  included  in  section  2  and  summarized  in  section  6.  The  literature  review  is  not 
intended  to  be  exhaustive.  This  report  highlights  the  main  issues  and  trends  in  the  fatigue 
research,  which  are  potentially  useful  for  improved  understanding  of  the  LCF  performance  of 
Army  land  combat  structures. 


2.  LCF-Related  Issues  in  Army  Applications 

Many  composite  structures  used  in  Army  applications  are  subjected  to  severe  fatigue  loading 
conditions  and  harsh  environment  over  long  periods  of  time.  Their  fatigue  life  and  performance 
depend  on  their  fatigue  damage  tolerance  and  the  rates  of  property  degradation  under  cyclic  or 
pulse  loads.  The  LCF  conditions  occur  when  an  Army  land  combat  system  is  used  under  high 
loads  that  may  be  as  high  as  90-95%  of  the  ultimate  strength  (section  2.2).  Such  high  values  of 
loading  are  reached  in  a  fraction  of  a  second.  The  low  cycle  and  high  cycle  (i.e.,  lower  load  level) 
fatigue  have  common  damage  mechanisms  and  similar  property  degradation  mechanisms.  Both 
LCF  and  high  cycle  fatigue  (HCF)  conditions  are  relevant  to  the  following  examples  of  Army 
applications: 

•  composite  armored  vehicles  (e.g.,  Crusader,  Scout,  vehicle  shock  absorption  systems,  and 
next  generation  vehicles  [Davila,  Chen,  and  Baker  1998]); 

•  composite  assault  bridge  and  future  line-of-communication  bridging; 

•  weapon  systems  (e.g.,  Howitzer  trails,  mortar  base  plate,  electromagnetic  gun  components, 
Crusader  cannon,  gun  tube  overwraps,  etc.  [Tzeng  1999]); 

•  rotorcraft  systems  (e.g.,  Comanche,  Apache,  etc.  [O'Brien  et  al.  1998]); 
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•  lightweight  personnel  protection  systems. 

In  the  following  section,  the  material  issues  relevant  to  fatigue  performance  of  composite 
structures  used  in  various  Army  land  combat  systems  are  discussed.  First,  examples  of  design 
requirements  for  fatigue  performance  in  several  Army  applications  are  discussed  in  section  2.1. 
A  brief  review  of  existing  fatigue  analysis  codes  and  complementary  analysis  tools,  which  can  be 
used  for  structural  and  microstructural  design  optimization,  is  also  presented  in  section  2.1. 
Most  common  industrial  applications  of  LCF  models  are  highlighted  in  section  2.2.  Fatigue 
issues  in  Army  applications  involving  the  thick-section  composite  are  analyzed  in  section  2.3. 
The  role  of  interphases  and  adhesive  joints  on  fatigue  performance  is  discussed  in  sections  2.4 
and  2.5,  respectively.  Environmental  effects  on  Army  land  combat  systems  are  reviewed  in 
section  2.6.  Recommendations  for  fundamental  research  in  the  areas  critical  for  understanding  the 
LCF  phenomena  are  included. 

2.1  Design  for  Fatigue  Performance.  The  fatigue  performance  characteristics  of  the 
composite  materials  used  are  required  for  finding  the  optimal  design  solutions  for  numerous 
Army  structures.  The  following  sections  highlight  software  solutions  that  can  be  used  to  meet 
the  design  and  fatigue  analysis  needs.  For  design  of  Army  systems,  special  attention  is  given  to 
the  effects  of  high  loads  on  the  thick-section  composite  structures  in  the  LCF  conditions. 
Postrepair  performance  of  new  hybrid  composites  is  also  evaluated  to  assess  the  retained  ballistic 
protection  capabilities. 

Knowledge  of  the  property  degradation  rates  is  especially  important  for  the  design  of 
systems  in  which  fatigue  performance  is  critical.  The  design  requirements  for  fatigue 
performance  vary  from  one  Army  system  to  another.  For  instance, 

•  the  composite  base  plate  of  a  mortar  should  withstand  the  maximum  pressure  of  about  16  ksi 
for  at  least  3,000  service  cycles, 

•  the  composite  barrel  of  a  mortar  is  designed  for  a  similar  range  of  pressures  but  for  104  fatigue 
loading  cycles, 

•  the  composite  compulsator  for  an  electromagnetic  gun  is  design  to  be  in  service  for  about  104 
loading  cycles  (in  service,  the  rotating  gun  compulsator  experiences  strain  rates  between  1  and 
10  (in/in  s)  while  it  decelerates;  the  peak  banding  stress  is  calculated  to  be  155  ksi  for  the 
maximum  rotor  tip  speed  of 400  m/s),  and 

•  the  howitzer  trails  are  designed  for  5,000  effective  full  charges  under  maximum  loads. 
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To  address  the  Army  after  next  (AAN)  goals,  most  of  the  aforementioned  design  requirements 
should  be  doubled  or  the  weight  of  structures  should  be  halved. 

Most  of  the  design  decisions  are  currently  based  on  the  static  mechanical  properties  of 
composite  materials  and  simple  relations  for  fatigue  performance.  Such  phenomenological  criteria 
are  predominantly  based  on  experimental  data  and  have  limited  predictive  capabilities. 
Limitations  of  current  fatigue  life  prediction  models  for  composites  force  large  factors  of  safety 
for  the  designed  structures.  The  associated  deficiencies  in  the  design-for-fatigue  methodologies 
lead  to  heavier  “overdesigned”  structures  that  tend  to  be  more  costly  than  necessary. 

2.1.1  Fatigue  Simulation  Codes  for  Design.  The  MRLife  simulation  code  (Case  and 
Reifsnider  1998)  is  developed  for  the  simulation  of  performance  and  fatigue  life  prediction  of 
composite  laminates.  MRLife  is  suited  for  fatigue  analysis  of  critical  material  elements  in  a  wide 
variety  of  problems  involving  polymer  matrix  systems  (see  Appendix).  Such  problems  include 
delamination  and  failure  of  notched  and  unnotched  materials  with  or  without  moisture  diffusion. 
The  pointwise  effects  caused  by  thermal  loads,  creep,  stress  relaxation,  and  aging  can  also  be 
accounted  for  in  the  MRLife's  fatigue  analysis  of  material  elements.  The  LCF  behavior  of 
polymer  matrix  composites  (PMCs)  can  be  also  predicted  by  using  the  MRLife  fatigue 
performance  simulation  code. 

Engineering  Mechanics  Technology,  Inc.  (Harris  and  Dedhia  1997)  has  developed  two  codes, 
Smart  Crack  and  Non-Linear  Smart  Crack,  for  the  linear  and  nonlinear  analysis  of  fatigue  crack 
propagation  in  metals.  The  linear  analysis  is  based  on  linear  elastic  fracture  mechanics  and  use  of 
the  mode  I  stress  intensity  factor  for  characterization  of  the  near-crack-tip  conditions.  In  the 
nonlinear  analysis,  the  plasticity  is  described  by  the  Ramberg-Osgood  constitutive  relation.  The 
creep  is  simulated  by  a  power  law  relation.  The  capability  of  the  metal-oriented  codes  to 
perform  fatigue  analysis  of  composite  structures  is  obviously  limited. 

The  impact  of  various  microstructures  on  macroscopic  behavior  of  composites  can  be 
analyzed  by  using  the  theory  of  homogenization  (Ericksen  1986)  or  “smearing-unsmearing” 
methodology  (Chou,  Carleone,  and  Hsu  1972).  This  methodology  is  implemented  in  a  software 
program  called  LAMP  AT,  which  was  developed  by  Bogetti,  Hoppel,  and  Bums  (1995),  for 
structural  analysis  and  design  of  thick-section  laminated  composite  structures.  The  code 
LAMP  AT  is  intended  to  be  used  in  conjunction  with  typical  commercial  finite  element  analysis 
(FEA)  codes.  At  the  preprocessing  stage,  LAMP  AT  is  utilized  to  generate  the  effective, 
homogeneous,  three-dimensional  properties  of  a  laminate,  which  are  used  as  input  for  an  FEA 
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code.  It  may  be  possible  to  link  the  LAMP  AT  capabilities  with  a  fatigue  performance  simulation 
code  (see  Appendix). 

For  the  structural  design  and  stress  analysis  purposes,  LAMP  AT  can  be  employed  as  a 
postprocessing  tool.  It  can  conduct  a  detailed  ply-level-based  failure  assessment  of  a  composite 
structure.  The  failure  assessment  is  based  on  a  wide  variety  of  lamina  failure  criteria.  The  results 
can  be  portrayed  graphically  by  PATRAN  (PDA  Engineering)  in  order  to  visualize  the  critically 
loaded  regions  within  a  structure.  A  postprocessing  program  with  similar  capabilities  has  been 
developed  by  Harik  (1997)  to  carry  out  complementary  analyses  of  stress  and  strain 
distributions  obtained  from  the  FEA  of  composite  structures  performed  by  the  FEA  commercial 
code  ABAQUS  (Hibbitt,  Karlsson  and  Sorensen  Inc).  The  program  is  called  Post-ABAQUS,  as 
it  is  designed  to  extend  the  postprocessing  capabilities  of  the  commercial  code  ABAQUS-POST. 
These  postprocessing  programs  can  be  supplemented  with  simple  S-N  relations  (i.e.,  stress  level 
vs.  number  of  cycles)  for  the  basic  fatigue  analysis. 

Another  post-FEA  program,  called  MSC/FATIGUE,  has  been  developed  by  MacNeal- 
Schwendler  Co.  It  offers  new  simulation  capabilities  for  the  effects  of  random  vibrations  on 
material  properties.  The  fatigue  life  predictions  are  based  on  frequency  response  and  the  random 
vibration  FEA  techniques.  The  program  is  especially  useful  for  the  design  of  electronic 
components,  microelectromechanical  systems  (MEMS),  wind  turbines,  and  various  components 
in  engines. 

2.1.2  Recommendations  for  Fatigue  Design.  Development  of  new  design  methodologies 
that  would  include  the  assessment  of  fatigue  performance  is  closely  linked  with  the  development 
of  new  fatigue  modeling  capabilities  and  associated  design  strategies.  The  short-term  design 
requirements  can  be  met  by  experimental  characterization  of  specific  composite  materials  (i.e.,  S- 
N  curves,  residual  strength  analysis,  stiffness  degradation  rates,  etc.).  The  experimental 
characterization  should  be  preceded  by  thorough  stress  analysis  of  the  composite  structure 
involved  in  order  to  determine  the  critical  structural  elements,  the  types  of  critical  loads,  and 
expected  damage  and  failure  modes. 

The  experimental  data  can  be  used  to  develop  simple  fatigue  life  prediction  models  for 
specific  composite  structures.  It  should  be  noted  that  such  models  are  phenomenological  in 
nature  and  are  limited  to  specific  groups  of  composite  materials  and  types  of  loading  under 
consideration  (see  section  5.2).  Phenomenological  fatigue  models  can  be  linked  with  LAMPAT's 
structural  analysis  capabilities  (Bogetti,  Hoppel,  and  Bums  1995).  The  corresponding  design- 
for-fatigue  strategies  may  also  include  other  postprocessing  software  tools. 
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In  order  to  meet  the  long-term  design  requirements  and  carry  out  multilevel  structural 
optimization  analyses,  one  needs  to  overcome  the  limitations  of  simple  phenomenological  models 
and  acquire  advanced  modeling  capabilities  for  fatigue  performance  assessment.  To  avoid 
excessive  duplication  on  one  hand  and  accelerate  the  progress  in  the  design-for-fatigue  efforts  on 
the  other  hand,  it  is  important  to  consider  and  evaluate  current  state-of-the-art  fatigue  simulation 
technologies.  The  MRLife  simulation  code  (see  Appendix)  is  suited  for  fatigue  performance 
analysis  of  a  wide  variety  of  polymer  and  ceramic  composite  systems.  The  MRLife  code  may 
simulate  delamination  and  failure  of  notched  and  unnotched  materials  with  or  without  moisture 
diffusion.  The  effects  caused  by  thermal  loads,  creep,  stress  relaxation,  and  aging  can  also  be 
accounted  for  in  MRLife's  fatigue  analysis.  The  LCF  modeling  capabilities  of  the  MRLife  code 
can  be  validated  via  selected  benchmark  experiments. 

The  modeling  capabilities  of  the  MRLife  fatigue  performance  simulation  code  can  be  linked 
with  the  structural  analysis  capabilities  of  commercial  FEA  codes  (e.g.,  ANSYS,  ABAQUS,  etc.). 
The  postprocessing  software  programs  like  LAMP  AT  can  be  incorporated  into  the  design-for- 
fatigue  strategies  for  complementary  stress  and  strain  analyses.  Development  of  more  robust 
methodologies  for  the  assessment  of  long-term  properties  should  result  in  better  design  strategies 
and  more  efficient  use  of  composite  materials  in  weight-critical  applications. 

2.2  LCF  Behavior  of  Composites.  LCF  is  a  critical  loading  condition  for  many  Army  land 
combat  composite  systems  and  is  not  well  understood.  The  unique  features  of  LCF  effects  on 
hybrid  composite  systems  have  not  been  studied  in  any  great  detail.  Most  of  publications  on  the 
LCF  involve  either  metals,  metal  alloys,  or,  in  some  cases,  metal  matrix  composites  (Coffin  1969; 
Solomon  et  al.  1987;  Rie  and  Portella  1998).  Corresponding  applications  include: 

•  gas  turbine  and  engine  components  made  of  titanium  alloys  for  high-temperature  creep-fatigue 
conditions  (e.g.,  compressor  discs,  the  disc  rim,  rotor  blades,  etc.), 

•  pressure  vessels  in  power  industries  (e.g.,  cast-iron  combustion  chambers,  steel  pipes  and 
tubes,  etc.), 

•  multicycle  forming  operations  (e.g.,  forging,  rolling  and  drawing), 

•  defense  applications  (helicopter  gears,  shock  absorbers,  etc.),  and 

•  structural  applications  in  construction  industry  (Barnes  1990). 
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The  LCF  of  PMCs  or  thick-section  hybrid  structures  has  received  little  attention.  Walrath 
and  Adams  (1980)  examined  the  LCF  range  for  the  pure  epoxy  resin  during  standard  fatigue 
characterization  tests  in  late  1970s.  Other  researchers  have  also  reached  the  LCF  loads  in  some 
of  their  fatigue  experiments  with  PMCs;  however,  the  unique  features  of  the  LCF  range  have  not 
been  identified  (Mandell  1981;  Barnes  1990;  Reifsnider  1991;  Case  and  Reifsnider  1998). 
Chaphalkar  (1998)  and  Harik  et  al.  (1999)  are  among  few  researchers  who  have  investigated  the 
effects  of  high  LCF  loads  on  thick-section  PMCs  in  work  recently  performed  under  the  Army- 
sponsored  Tuskegee  Research  Consortium  and  ARL  Postdoctoral  Research  Program. 

2.2.1  Basic  Physical  Phenomena  Relevant  to  the  LCF  Studies.  LCF  conditions  are 
characterized  by  the  high  loads  applied  to  a  material  over  a  fraction  of  a  second  and  then  repeated 
for  a  finite  number  of  loading  cycles.  The  corresponding  strain  rates  may  vary  between  1  and  10 
in/in-s,  or  even  more  than  100  in/in-s  for  some  weapon  systems.  Since  the  LCF  loads  may  reach 
90-95%  of  the  static  ultimate  strength,  the  material  behavior  of  the  composite's  constituents  may 
no  longer  be  regarded  as  linear.  Indeed,  when  the  strains  exceed  4-5%,  the  matrix  material  and  the 
interphase  undergo  finite  deformation,  while  fibers  may  experience  interfacial  debonding  and 
fracture.  The  strain  rates  and  nonlinear  deformation  encountered  under  a  LCF  loading  are  not  as 
high  as  those  seen  during  a  ballistic  impact.  However,  it  is  plausible  that  the  material  degradation 
mechanisms  and  material  properties  affecting  the  LCF  performance  of  composites  are  similar  to 
those  determining  the  ballistic  performance  of  Army  combat  composite  systems.  The  extent  of 
possible  correlation  between  the  LCF  and  ballistic  performance  of  composite  materials  is  not 
known. 

Section  3  presents  a  review  of  typical  fatigue  damage  mechanisms  that  are  activated  in 
composites  having  different  microstructures.  The  effect  of  LCF  conditions  on  typical  damage 
mechanisms  needs  to  be  thoroughly  examined.  Patterns  of  damage  evolution  in  the  laminates 
with  typical  structures  (e.g.,  unidirectional,  cross-ply,  and  woven  fiber  composites)  are  also 
discussed  in  section  3.  Classification  of  different  stages  in  the  damage  accumulation  processes  is 
presented  in  section  4.  A  summary  of  the  damage  accumulation  models  and  modeling  approaches 
is  included  in  section  5.  The  effects  of  LCF  conditions  on  damage  accumulation  processes  and 
the  property  degradation  rates  need  to  be  evaluated. 

2.2.2  Recommendations  for  LCF  Studies.  More  fundamental  research  is  needed  in  the  area  of 
modeling  and  characterization  of  polymer  composite  materials  subjected  to  low-cycle  or  pulse- 
vibration  fatigue  (PVF).  The  need  to  develop  a  physical  understanding  of  the  mechanisms  and 
phenomena  associated  with  fatigue-induced  failure  in  PMCs  is  especially  important.  The  effects 
of  LCF  conditions  on  typical  damage  mechanisms  and  the  extent  of  plastic  deformation  (Hill 
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1950;  Drucker  1967)  can  be  evaluated  by  thorough  experimental  testing  and  detailed 
nondestructive  monitoring  of  fatigue  damage  accumulation  processes  (Sasaki  1 997).  Comparative 
analysis  of  material  degradation  mechanisms  present  during  LCF  and  ballistic  loadings  is  required 
in  order  to  determine  the  extent  of  possible  correlation  between  the  LCF  and  ballistic 
performance  of  composite  structures  in  Army  land  combat  applications.  Experimental 
characterization  is  also  needed  for  the  development  of  novel  LCF/PVF  models  and  realistic 
service  life  prediction  of  thick-section  composites  in  various  Army  systems.  More  accurate 
service  life  prediction  capability  for  Army  structures  under  LCF  conditions  would  enable 
optimization  for  performance  and  cost. 

2.3  LCF  of  Thick-Section  Composites.  The  LCF  behavior  of  thick-section  composites  is 
of  special  interest  to  Army  applications.  Multi-ply  composites  are  part  of  integral  armor 
composite  systems.  Advances  in  integral  armor  and  other  complex  composite  systems  depend 
on  a  thorough  understanding  of  possible  failure  mechanisms.  Although,  mechanics  of  anisotropic 
plates  (Lekhnitskii  1968)  and  layered  materials  (Chou,  Carleone,  and  Hsu  1972;  Christensen 
1979;  Ericksen  1986)  has  been  vigorously  investigated,  general  fatigue  behavior  and  LCF  of  thick- 
section  composites  have  received  little  attention  to  date  (Solomon  and  Halford  1987;  Rie  and 
Portella  1998).  LCF  behavior  of  homogeneous  thick- section  composites  must  be  understood 
before  LCF  of  their  heterogeneous  hybrid  counterparts  can  be  addressed.  A  review  of  typical 
mechanisms  of  fatigue  damage  accumulation  is  presented  in  section  3.  In  particular,  the 
dependence  of  patterns  of  damage  evolution  on  different  microstructures  is  analyzed. 
Understanding  the  fatigue  damage  mechanisms  and  their  interaction  with  overall  mechanical 
response  is  needed  to  adequately  model  and  predict  LCF  behavior. 

2.3.1  LCF  Testing  of  Thick-Section  Composites.  LCF  performance  of  multilayer  composites 
has  not  been  studied  in  any  great  detail  either  analytically  or  experimentally.  Fatigue  behavior  of 
composite  materials  is  typically  evaluated  by  experiments  involving  either  tension-tension  or 
tension-compression  loading.  Even  for  non-hybrid  thick-section  composites,  sufficient  load 
transfer  between  tensile  grips  and  the  specimen  is  difficult  to  achieve  (Bakis  and  Stinchcomb 
1986;  Chaphalkar  1998).  During  flexural  fatigue,  which  is  relevant  to  the  thick-section  integral 
armor  applications,  both  compressive  and  tensile  loads  are  present.  Bending  fixtures  are 
sufficiently  robust  for  LCF  testing  of  thick-section  hybrid  composite  beams. 

Composite  materials  with  toughened  matrices  usually  have  better  fatigue  performance.  The 
mechanisms  of  toughening  in  PMCs  include  interfacial  debonding  and  interfacial  void  formation 
around  rubber  or  thermoplastic  particles,  plastic  shear  localization,  and  dissipation  of  energy 
through  viscoelastic  relaxation.  Viscoelastic  relaxation  is  accompanied  by  a  temperature  rise; 
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because  the  polymer  matrix  is  a  poor  thermal  conductor,  the  dissipation  of  heat  is  slow  and  local 
temperature  increases  (Joseph  1990).  Shen,  Chen  and  Sauer  (1983)  demonstrated  that  at  higher 
stresses  hysteretic  heating  is  quite  extensive  in  the  specimens  of  cast  poly-methyl  methacrylate 
subjected  to  tension-compression  fatigue.  As  a  result,  so-called  hot  spots,  extensive  plastic 
deformation  and  voids  are  possible  in  PMCs.  Thus,  premature  catastrophic  failure  of  a 
composite  structure  can  result. 

Fundamental  research  is  required  to  determine  the  role  of  localized  heating  and  viscoelastic 
interfaces  (Holmes  et  al.  1999)  in  the  failure  of  composites  under  LCF  conditions.  Development 
of  robust  numerical  techniques  for  viscoelastic  models  with  geometric  nonlinearities  is  especially 
important  (Harik  1997;  Harik  and  Caimcross  1999).  A  better  understanding  of  the  viscoelastic 
failure  mechanisms  will  allow  for  the  improved  ability  to  predict  catastrophic  failures  of 
composites  under  fatigue  through  computer  modeling. 

2.3.2  Recommendations  for  Studies  of  Thick-Section  Composites.  The  load  transfer 
problems  commonly  encountered  during  tension  and  compression  testing  of  thick-section 
composites  are  easily  avoided  in  flexural  fatigue  testing.  Flexural  fatigue  tests  can  be  carried  out 
by  using  a  three-point  or  a  four-point  bending  fixture  scaled  up  from  existing  testing  fixtures 
(Whitney,  Daniel,  and  Pipes  1984;  Carlsson  and  Gillespie  1990).  Recent  advances  in 
nonintrusive  techniques  for  nondestructive  testing  allow  systematic  analysis  of  fatigue  damage 
accumulation  within  the  specimen  being  tested.  Such  techniques  include  x-rays,  x-ray-based 
computed  tomography,  video-imaging,  acoustic  wave  scattering,  infrared  thermal  imaging,  and 
embedded  fiber  optics  (Carlsson  and  Gillespie  1990;  Fink  and  Corona-Bittick  1999;  Flores  et  al. 
1998).  Material  damage  assessment  prior  to,  during,  and  following  fatigue  tests  can  enable  greater 
understanding  of  material's  fatigue  behavior  and  damage  accumulation  processes.  From  this 
information,  possible  deviations  from  typical  damage  mechanisms  can  be  evaluated  and 
understanding  of  fatigue  performance  of  the  thick-section  composites  can  be  advanced. 

Comprehensive  analysis  of  the  fatigue  damage  data  (e.g.,  damage  location,  type,  shape,  size, 
and  evolution  patterns  [Talreja  1993;  Sasaki  1997])  and  gross  experimental  measurements  of 
macroscopic  mechanical  properties  (e.g.,  stiffness,  residual  strength)  can  provide  the  physical 
foundation  for  future  analytical  and  finite  element  fatigue  models.  Once  the  fatigue  behavior  of 
homogeneous  multi-ply  composites  is  more  fully  understood,  the  fatigue  of  heterogeneous  hybrid 
composites  can  be  analyzed  by  more  sophisticated  fatigue  models.  It  should  be  noted  that  the 
hybrid  composites  possess  layers  with  significant  differences  in  stiffnesses  that  would  lead  to 
the  deformations  violating  the  basic  assumptions  used  in  classical  bending  analyses. 
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Ultimately,  a  more  complete  understanding  of  fatigue  damage  mechanisms  in  thick-section 
and  integrated  composites  will  be  invaluable  to  advancing  the  applicability  of  composites  in  the 
Army  systems.  The  resulting  knowledge  incorporated  into  useful  models  should  yield  reliable 
predictions  of  the  fatigue  life  of  composite  structural  components  subjected  to  various  fatigue 
conditions.  It  will  enable  life-cycle  design  of  thick-section  composite  Army  structures  subjected 
to  LCF. 

2.4  Optimization  of  Interphases  for  LCF  Performance.  A  primary  concern  for 
composites  being  used  under  LCF  conditions,  which  involve  high  loads,  is  premature  failure. 
Thick-section  composites  typically  fail  at  stresses  and  strains  that  are  well  below  the  expected 
failure  limits.  This  early  failure  is  often  attributed  to  the  existence  of  critically  sized  processing 
and/or  material  defects  and  interfacial  problems  in  the  interphase  region  between  the  matrix  and 
the  reinforcing  phase  (Drzal  1983,  1986;  Sottos  1990;  Palmese  1992;  Skourlis  1995;  Hrivnak 
1996;  Harik  1997;  VanLandingham  1997;  Fink  and  McCullough  1999).  Evaluation  of  interphasial 
mechanical  properties  can  be  carried  out  experimentally  (Sottos  1990;  VanLandingham  1997)  or 
theoretically  (Palmese  1992;  Chu  and  Rokhlin  1996).  Modification  of  either  the  polymer  matrix 
or  the  reinforcing  fibers  to  improve  the  adhesion  between  the  material  components  has  proven  to 
be  a  key  in  the  optimization  of  performance  of  various  composite  structures.  Such  optimization 
should  be  based  on  rapidly  developing  interfacial  mechanics  of  fiber-reinforced  materials  (Clyne 
and  Watson  1991 ;  Harik  1997;  Harik  and  Lambros,  to  be  published). 

The  extent  of  the  interphase  region  in  composites  is  significant  (Hughes  1991).  For  instance, 
a  1  cm3  of  a  composite  is  filled  with  a  fiber  volume  content  of  60%  and  contains  as  many  as  3 
million  single  filaments.  The  total  area  of  the  fiber  surface  is  3,400  cm2.  As  a  result,  the  matrix 
and  its  ability  to  adhere  to  a  fiber  are  paramount  to  the  effective  transfer  of  the  mechanical  load  in 
the  composite  (Eriksonand  Plueddemann  1974;  Drzal  1983,  1986;  Fishman  1991;  Piggott  1991). 
The  interface  between  matrix  and  fiber  has  many  commonalties  with  the  interfacial  region  in 
laminated  systems  (O'Brien  1991).  In  both  cases,  a  large  surface  area  plays  a  direct  role  in  the 
load  transfer  from  the  matrix  to  the  reinforcing  constituent.  Stinchcomb  and  Reifsnider  (1975) 
emphasized  that  “the  way  the  interface  interacts  with  the  matrix  and  with  the  fibers  is  quite 
important  in  determining  fatigue  damage  initiation  in  composite  materials.”  Reifsnider  (1994), 
Subramanian,  Reifsnider,  and  Stinchcomb  (1995),  and  Subramanian  et  al.  (1996)  demonstrated 
that  the  fiber-matrix  interphase  can  significantly  affect  the  mechanical  properties  and  fatigue 
behavior  of  composites. 

The  fiber-matrix  interphases  (Drzal  1983;  Hughes  1991;  Palmese  1992;  Sottos,  Hiemstra,  and 
Scott  1994;  Hrivnak  1996;  Fink  and  McCullough  1999)  are  known  to  affect  the  local  material 
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properties  (Skourlis  1995;  VanLandingham  1997),  stress  distribution  (Sottos  1990;  Cervenka 
1995;  Kharik  1997),  interphasial  deformation  (Grabovsky  and  Kohn  1995;  Harik  and  Lambros, 
to  be  published)  and  nucleation  of  interfacial  cracks  (Harik  1997),  which  may  interact  with 
transverse  cracks  (Bailey  and  Parvizi  1981)  and  affect  the  fatigue  behavior  (Reifsnider  1994; 
Lesko,  Rau,  and  Riffle  1995).  Reinforcing  fibers  and  particles  themselves  may  serve  as  stress 
raisers  and  lead  to  interfacial  cracking  (Eshelby  1957).  Fiber-matrix  debonding  and  cracks  may 
significantly  reduce  the  load  transfer  between  matrix  and  the  fibers  and  cause  cracking  in 
composites  (Sottos,  Li,  and  Agrawal  1994;  Budiansky,  Hutchinson,  and  Slutsky  1995). 
Interfacial  damage  (Keer,  Dundurs,  and  Kiattikomol  1973;  Hashin  1991;  Pan,  Green,  and  Heilman 
1996)  or  material  inhomogeneity  of  interphases  also  affects  the  elastic  properties  of  composites 
(Jasiuk  and  Kouider  1993;  Lagache  et  al.  1994;  Low  et  al.  1995;  Theocaris  and  Demakos  1995; 
Lutz  and  Zimmerman  1996),  the  residual  stresses  (Jayaraman  and  Reifsnider  1993),  and  their 
macroscopic  behavior  (Tsai,  Arocho,  and  Gause  1990;  Kharik  1997;  Kim  and  Mai  1991, 1998). 

2.4.1  Fracture  Toughness  of  Fiber/Matrix  Interfaces.  In  understanding  failure  at  the 
interface,  one  must  closely  examine  the  polymer  matrix  and  its  interaction  with  the  interfacial 
surfaces  (Wool  1995;  Hrivnak  1996).  It  has  been  known  that  surface  treatments  can  improve  the 
interfacial  bond  strength  in  fiber-reinforced  composites  (Shorthall  and  Yip  1976).  Recently, 
research  on  polymer  films  demonstrated  that  a  polymer  chain  adapts  to  the  presence  of  an 
interface  so  the  physical  properties  of  the  polymer  films  are  affected  (Drzal  1986;  Palmese  and 
McCullough  1994;  Skourlis  1995;  Fink  and  McCullough  1999;  VanLandingham  et  al.  1999). 
Such  phenomena  may  lead  to  reduced  glass  transition  temperatures,  increased  diffusion  rates,  and 
nonuniform  curing  (Palmese  1992;  Skourlis  1995).  In  a  thermoplastic  polymer,  such  as 
polystyrene,  the  surface  has  a  glass  transition  temperature  that  is  significantly  lower  than  that  for 
the  bulk  material  (Skourlis  1995). 

Lesko,  Rau,  and  Riffle  (1995)  investigated  fatigue  performance  of  a  woven  carbon/vinyl  ester 
composite  that  had  carbon  fibers  sized  with  the  thermoplastic  coating.  The  sizing  considered  led 
to  an  improved  bond  between  the  matrix  and  fibers.  This  resulted  in  better  fatigue  durability  of 
the  composite  with  thermoplastic  interphase  in  relation  to  the  composites  with  unsized  carbon 
fibers.  Oyama  et  al.  (1996)  show  that  the  interdiffusion  at  the  interface  between  the  poly- vinyl 
pyrollidone  sizing  considered  by  Lesko,  Rau,  and  Riffle  (1995)  and  the  vinyl  ester  resin  is  easily 
facilitated.  This  leads  to  an  order  of  magnitude  larger  interphasial  thickness,  which  allows  greater 
energy  absorption  by  the  interphase. 

2.4.2  Recommendations  for  Studies  of  Interphases.  The  LCF  performance  of  composite 
structures  can  be  affected  by  improving  energy-absorbing  capabilities  of  the  fiber/matrix 
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interphase  (DiAnselmo,  Accorsi,  DiBenedetto  1992;  Kim  and  Mai  1998).  Subramanian  et  al. 
(1995,  1996)  demonstrated  that  the  interphasial  properties  may  influence  the  fatigue  behavior  of 
composites.  The  energy  absorption  mechanisms  such  as  fiber-matrix  debonding,  interfacial  void 
growth,  fiber  push-out,  or  frictional  fiber  sliding  (Kim  and  Mai  1998)  should  be  extensively 
investigated.  A  new  test  apparatus  called  Dynamic  Interphase-Loading  Apparatus  (DILA), 
developed  by  Tanoglu  et  al.  (to  be  published),  can  be  employed  to  characterize  the  fiber-matrix 
interphase  properties  around  an  individual  fiber  under  various  strain  rates.  This  knowledge 
should  also  contribute  to  the  understanding  of  the  ballistic  performance  of  composite  structures, 
develop  new  methodology  for  experimental  characterization  of  interphase  mechanical  properties 
(Wang  and  Chiang  1996),  and  advance  the  physical  understanding  of  various  interfacial 
phenomena  (Verpoest  and  Jones  1991). 

The  AAN  armor  systems  may  benefit  from  the  hybrid  composite  structures  involving  glass 
and  carbon  fibers.  Such  composite  systems  would  have  high  durability  and  corrosion  resistance 
due  to  glass  fibers  and  superior  specific  stiffness  and  strength  due  to  the  carbon  fiber 
reinforcement  (Mahfuz  et  al.  1995,  1996).  It  is  known  that  glass  fiber/vinyl  ester  composites 
have  good  static  and  long-term  fatigue  durability.  The  majority  of  commercially  available  carbon 
fiber  tows  are  not  sized  for  optimal  chemical  compatibility  with  the  vinyl  ester  resin.  Therefore, 
the  adequacy  of  the  strength  and  fatigue  performance  of  composites  with  the  carbon  fiber/vinyl- 
ester  interphases  is  uncertain.  The  fatigue  characterization  of  vinyl-ester  systems  is  also 
important  because  of  advances  in  new  co-injection  resin  transfer  molding  (CIRTM)  technology 
involving  vinyl-ester  resin  (Fink,  McKnight,  and  Gillespie  1998;  Fink  and  Gillespie  1999). 

2.5  Fatigue  of  Adhesive  Joints  in  Composite  Structures.  Composite  integral  armor  can 
be  viewed  simplistically  as  a  composite  system  of  five  basic  subsystems  that  are  connected  by 
adhesive  bonds.  The  composite  subsystems  include  a  thin  graphite/epoxy  or  glass/epoxy  face 
sheet,  a  ceramic  tile  for  ballistic  protection,  a  rubber  layer  to  mitigate  damage  to  the  composite 
backing  during  initial  fracture  of  the  ceramic,  and  a  thick  composite  layer  for  structural  and 
ballistic  performance.  Each  layer  has  very  distinctive  mechanical  properties  that  induce  unique 
stress  distributions  in  the  adhesive  layers.  The  prevalent  loading  modes  are  compression  and 
flexure  (Davila,  Chen,  and  Baker  1998).  The  integrity  of  adhesive  joints  and  interfaces  must  be 
maintained  under  severe  fatigue  loading  in  order  to  retain  the  ballistic  performance.  It  is 
important  before  and  after  a  ballistic  impact  and  after  consequent  repairs.  Appropriate  analysis 
can  be  carried  out  with  analytical  models  that  take  into  account  the  effect  of  inclusions  on  the 
mechanical  response  of  plates  subjected  to  poit  loads. 
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The  adhesion  of  the  interphases  between  the  constituent  materials  strongly  affects  ballistic 
performance,  damage  tolerance,  and  long-term  durability.  Modification  of  these  interfaces  could 
change  the  modes  of  failure  and  energy-absorption  characteristics.  To  optimize  the  ballistic 
performance,  it  is  important  to  isolate  the  effect  of  high-strain-rate  loading  in  the  adhesive 
bondline.  The  interfaces  between  dissimilar  materials  in  hybrid  composite  structures  (Delale  and 
Erdogan  1988)  have  many  common  features  with  the  interphase  between  the  matrix  and  a  fiber. 
Indeed,  both  cases  involve  large  surface  areas  that  play  an  important  role  in  the  load  transfer  from 
one  constituent  material  to  another.  Henceforth,  the  fatigue  performance  studies  of  both  types  of 
interphases  can  be  linked  together.  New  processing  techniques,  such  as  CIRTM  (Fink, 
McKnight,  and  Gillespie  1998;  Fink  and  Gillespie  1999)  that  has  been  shown  to  provide 
improved  ballistic  properties  in  composite  armor  and  create  diffused  interphases  between 
dissimilar  resin  layers  in  a  composite  structure  (Fink,  McKnight,  and  Gillespie  1998). 

2.5.1  Repair  of  Composite  Structures.  A  ballistic  impact  of  a  composite  armor  structure 
induces  multiple  cracks  and  delamination  within  the  structure.  Interfacial  debonding  and 
delamination  fractures  are  effective  mechanisms  for  absorption  of  projectile  energy  (Monib  et  al. 
1999).  However,  the  defects  induced  may  reduce  the  residual  strength  of  the  composite  structure 
so  that  it  may  no  longer  possess  the  required  structural  and  ballistic  characteristics.  If  such 
characteristics  are  below  the  design  allowables  then  immediate  repairs  are  needed.  Fundamental 
studies  are  needed  to  predict  the  postrepair  performance  of  hybrid  composite  structures  (Fink 
and  Gillespie  1999).  The  structural  integrity  of  repaired  components  may  influence  the 
constraining  pressure  effects  on  adjacent  regions  of  a  composite  structure  and  their  property 
degradation  (Coffin  and  Rogers  1967;  Harik  1997;  Harik  and  Caimcross  1999). 

The  repair  of  complex  composite  structures  should  be  carried  out  as  a  comprehensive 
structural  maintenance  program.  Establishing  such  a  program  should  include  evaluation  of  the 
following  aspects: 

•  sources  of  structural  deterioration  in  Army  systems, 

•  susceptibility  of  each  Army  structure  to  various  sources  of  deterioration, 

•  the  consequences  of  deterioration  to  continued  battle  worthiness, 

•  the  effectiveness  of  detection  methods  in  finding  structural  deterioration  (taking  into 
account  inspection  thresholds  and  intervals, 
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•  the  effectiveness  of  the  repair  in  restoring  load-carrying  and  ballistic  capabilities  and  the 
effect  on  the  integrity  of  an  Army  structure,  and 

•  the  effectiveness  of  prevention  and  control  measures  to  mitigate  existing  and  anticipated 
problems. 

2.5.2  Recommendations  for  Studies  of  Adhesive  Joints.  A  comprehensive  study  of  primary 
LCF  damage  evolution  mechanisms  in  ballistically  impacted  and  repaired  Army  structures  is 
needed.  Reliable  prediction  of  the  ballistic  damage  tolerance  of  various  Army  systems  is 
possible,  only  if  the  LCF  mechanics  of  damaged  hybrid  structures  is  well  understood.  The  LCF 
performance  of  damaged  composites  has  to  be  evaluated  in  order  to  set  the  guidelines  for  the 
repair  program.  The  postrepair  mechanics  of  complex  composite'systems  is  also  closely  linked 
with  the  mechanics  of  in-service  degradation  processes  such  as  fatigue,  corrosion,  and  the  effects 
of  discrete  ballistic  impacts. 

2.6  Environmental  Effects  on  Army  Systems.  Future  Army  applications  of  polymer- 
matrix  and  hybrid  composites  presume  that  the  composite  structures  will  experience  considerable 
exposure  to  a  wide  range  of  temperatures  and  humidity  levels.  Hence,  during  service  life,  the 
mechanical  fatigue  loading  is  often  complemented  with  various  environmental  conditions.  The 
resulting  physical  environment  may  involve  thermomechanical  loads,  moisture  effects,  chemical 
corrosion,  mechanical  wear,  etc.  In  the  moist  or  chemically  active  environment,  degradation  of 
mechanical  properties  may  be  especially  severe  since  the  pre-existing  defects  and  developing 
matrix  cracking  may  lead  to  significant  plasticization  of  the  matrix.  Dramatic  changes  in  the 
material  damage  state  can  lead  to  premature  failure. 

Adams,  Bowles,  and  Herakovich  (1986)  demonstrated  that  transverse  tracking,  which  is 
common  under  fatigue  conditions,  significantly  reduces  the  coefficient  of  thermal  expansion.  In 
the  1980s,  a  number  of  researchers*  showed  how  the  moisture  uptake  increases  due  to  matrix 
cracking  and  affects  the  matrix  material  properties.  Wang,  Bogetti,  and  Gillespie  (1998)  showed 
that  the  fiber-matrix  interphase  has  higher  moisture  content  that  leads  to  increased  stress 
concentration  and  higher  probability  of  interfacial  failure. 

2.6.1  Environmental  Effects  on  Integrity  of  Interfaces.  Interfacial  mechanical  properties 
may  deteriorate  due  to  a  number  of  environmental  effects.  In  humid  environments,  the  composite 


’  See  Wang,  Bogetti,  and  Gillespie  (1998). 
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structures  can  absorb  water  or  chemically  active  vapors.  Absorbed  moisture  can  affect  the 
residual  stress  states,  stress  concentration,  material  damage  state,  molecular  network  structure, 
and  glass  transition  temperature.  Furthermore,  plasticization  of  the  interphase  may  also  be 
caused  by  moisture  absorption  (Wang,  Bogetti,  and  Gillespie  1998).  Such  interphasial  changes 
would  affect  the  load  transfer  between  the  matrix  and  reinforcing  fibers  and  the  macroscopic 
damage  state.  In  the  toughened  matrices,  along  with  the  aforementioned  interfacial  phenomena, 
there  may  be  additional  interfacial  problems  such  as  crazing  around  rubber  particles. 

Fiber-reinforced  plastics  exhibit  values  of  hardness  and  thermal  resistance  that  are  lower  then 
those  of  metallic  materials.  Surface  microcracks  may  nucleate  from  the  surface  flaws  and  notches 
and  evolve  into  matrix  cracks  under  mechanical  and  environmental  fatigue  conditions.  As  the 
microcracks  grow,  they  may  interact  with  the  matrix  macrocracks.  This  would  increase  the  crack 
density  and  the  probability  of  changes  in  the  material  damage  state.  Protective  coating  of 
composites  (Paesano,  Visconti,  and  Penasa  1992)  may  impede  the  growth  of  surface  defects. 
Metallic  or  ceramic  coatings  usually  improve  the  hardness  and  thermal  fatigue  resistance  of 
PMCs  (Vishwanath,  Varna,  and  Rao  1990).  Such  coatings  have  advanced  surface  properties  and 
a  coat-substrate  adhesion  that  is  strong  enough  to  withstand  high  loads  present  in  structures 
under  LCF.  The  only  drawback  is  that  high  adhesion  is  achieved  by  high-temperature  processing, 
that  is  not  always  well  tolerated  by  the  plastic  substrate. 

2.6.2  Recommendations  for  Studies  of  Environmental  Effects.  Long-term  exposure  of  Army 
composite  structures  to  a  humid  and/or  chemically  active  environment  can  significantly  affect  the 
material  properties,  material  damage  states,  ballistic  performance,  and  fatigue  life  of  Army 
systems.  During  moisture  absorption,  various  fatigue  events  such  as  matrix  microcracking, 
interfacial  debonding,  and  property  degradation  of  the  matrix  and  the  interphase  take  place.  On 
the  other  hand,  the  external  mechanical  loads  may  initiate  different  structural  defects  and  induce 
the  material  damage  states  that  cause  higher  rates  of  moisture  absorption.  These  deleterious 
effects  require  in-depth  studies  especially  in  the  area  of  interfacial  properties  of  composites, 
kinetics  of  fatigue  damage,  and  material  damage  states.  The  microstructural  and  macrostructural 
interphases  play  a  critical  role  in  load  transfer  between  dissimilar  constituent  materials  and 
structural  components. 


3.  Fatigue  Damage  Mechanisms:  Structural  Effects 

Advances  in  fatigue  life  prediction  capability  require  understanding  of  the  fundamental 
mechanisms  of  material  degradation  during  service  life.  The  changes  in  macroscopic  mechanical 
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behavior  of  composites  and  the  corresponding  rates  of  property  degradation  strongly  depend  on 
the  patterns  of  damage  evolution.  High  loads  and  moderate  strain  rates  associated  with  the  LCF 
conditions  may  affect  the  initiation  of  damage  and  the  kinetics  of  fatigue  damage  accumulation 
between  different  material  damage  states.  Under  fatigue  loading  conditions,  composite  systems 
accumulate  extensive  microstructural  and  macroscopic  damage.  The  complex  nature  of  damage  in 
composite  systems  can  be  understood  only  if  the  microstructural  effects  are  also  considered,  as 
different  microstructures  lead  to  different  patterns  of  fracture  and  damage  accumulation  (Talreja 
1993, 1994),  which  result  in  different  material  damage  states  (Case  and  Reifsnider  1998). 

On  a  microscopic  level,  the  mechanical  damage  occurs  through  various  uncoupled  (i.e., 
noninteractive)  and  coupled  (i.e.,  interactive)  damage  and  failure  modes.  A  particular  damage 
mechanism  may  involve  a  number  of  different  damage  modes  (Talreja  1987).  On  a  macroscopic 
level,  accumulation  of  damage  leads  to  significant  changes  in  material  properties  (e.g.,  stiffness, 
strength,  thermal,  and  electric  conductivities  [Case  and  Reifsnider  1998]),  as  well  as  changes  in 
the  material  damage  state.  The  effects  of  high  loads  and  finite  strain  rates  on  various  damage  and 
failure  modes  depend  on  microstructural  elements  (e.g.,  properties  of  fibers  and  their  sizing, 
properties,  and  structure  of  interphases)  and  macrostructural  elements  (e.g.,  stacking  sequence, 
interlaminar  layers,  and  through-thickness  reinforcement). 

Laminated  composite  structures  have  inherently  low  interlaminar  fracture  toughness 
(Reifsnider  1991).  As  a  result,  their  most  prevalent  fatigue  failure  mode  is  delamination.  The 
process  of  interlaminar  delamination  of  composites,  its  initiation,  and  growth  are  complex 
mechanical  phenomena  (O'Brien  et  al.  1982;  O'Brien  1991).  The  delamination  failure  is  usually 
preceded  by  extensive  damage  accumulation  in  the  matrix,  fiber-matrix  interphase  and  interlaminar 
region  (Talreja  1987).  Different  types  of  fatigue  damage  influence  each  other  and  affect  the 
damage  evolution  patterns  (Broutmanand  Sahu  1969;  Mandell  1981;  Mandelll  et  al.  1984;  Pook 
1995, 1997, 1998;  Talreja,  to  be  published). 

The  evolution  of  fatigue  damage  in  most  composite  systems  can  be  described  by  a  series  of 
material  damage  states  (Reifsnider  1994)  and  a  set  of  damage  mechanisms  (Talreja  1987).  These 
mechanisms  can  be  divided  into  six  groups: 

•  evolution  of  existing  damage  (deformation  of  microvoids,  propagation  of  microcracks, 
void-crack  interaction,  coalescence  of  voids), 

•  fiber  damage  mechanisms  (evolution  of  surface  defects,  nucleation  of  surface  cracks  due  to 
surface  roughness,  failure  caused  by  the  interfacial  stress  concentration), 
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•  interfacial  damage  mechanisms  (nucleation  of  interfacial  cracks,  debonding  of  the  interface, 
interfacial  separation  and  sliding,  microcracking  in  the  interphase,  interfacial  void 
formation), 

•  matrix  damage  mechanisms  (nucleation  of  microcracks  from  existing  flaws,  propagation  of 
microcracks,  formation  of  transverse  macrocracks), 

•  interlaminar  damage  mechanisms  (interlaminar  debonding,  interlaminar  cracking),  and 

•  through-the-thickness  splitting  of  composite  laminates. 

As  noted  before,  a  particular  type  of  damage  mechanism  can  be  associated  with  one  or  more 
fracture  modes.  For  instance,  formation  of  a  transverse  matrix  microcrack  is  associated  with 
mode  I  crack  opening  that  is  caused  by  the  normal  tensile  stresses.  The  interfacial  crack  opening, 
which  may  affect  the  transverse  matrix  macrocracking,  is  associated  with  mode  I  and  mode  II 
crack  opening.  Mode  II  is  caused  by  the  pure  shear  failure.  Mode  I  has  dominant  influence  on 
the  interfacial  separation  (Hashin  1991;  Harik  1997).  Intrinsic  strength  of  a  composite  material  is 
often  controlled  by  these  highly  localized  stress  conditions  and  fracture  events. 

3.1  Fatigue  Damage  in  Unidirectional  Laminated  Composites.  The  activation  of 
particular  fatigue  damage  mechanisms  in  unidirectional  composites  depends  on  the  type  of 
loading  (e.g.,  tension-tension,  tension-compression,  compression-compression),  its  value  with 
respect  to  the  ultimate  strength  of  structural  components,  the  loading  rate,  and  the  direction  of 
loading  (e.g.,  parallel,  inclined  or  perpendicular  to  the  fiber  direction).  First,  the  fatigue  damage  is 
analyzed  under  the  most  common  tensile  loads  parallel  to  fibers. 

3.1.1  Damage  Caused  by  Axial  Fatigue  Loading.  Fatigue  damage  in  unidirectional 
composites  usually  starts  with  transverse  microcracking  in  the  matrix,  breakage  of  individual 
weak  fibers,  and  longitudinal  interfacial  debonding.  The  matrix  microcracking  can  evolve  into  a 
network  of  macrocracks  via  crack  propagation  and  crack  bridging.  The  evolution  of  macrocracks 
may  be  affected  by  the  interfacial  debonding  and  finite  separation  of  the  fiber-matrix  interface 
through  crack  interaction  and  crack  bridging.  The  surface  cracks  and  notches  may  evolve  into 
matrix  cracks  and  interact  with  the  matrix  macrocracks  as  well.  Initially,  matrix  cracking  can  be 
considered  homogeneous  and  non-interactive  throughout  individual  plies  (Talreja  1987).  The 
crack  density  may  reach  a  steady  state  as  the  characteristic  damage  state  (CDS)  sets  in 
(Reifsnider  1977).  The  transverse  cracks  may  also  induce  an  interlaminar  fracture  and  consequent 
delamination  of  different  plies  (see  sections  4.2  and  4.3).  This  stage  in  the  fatigue  damage 
accumulation  often  leads  to  the  final  failure  of  a  laminate  composite. 
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The  broken  fibers  cause  shear-stress  concentrations  at  the  interface  around  the  fracture  site 
(Gamstedt  and  Talreja  1999).  This  usually  leads  to  localized  interfacial  failure  with  the  length  of 
debonding  being  equal  to  a  few  fiber  diameters.  The  extent  of  debonding  depends  on  the  shear 
strength  of  the  interface.  The  segments  of  a  broken  fiber  usually  separate  and  form  a  void.  The 
void  formed  induces  stress  concentration  for  the  longitudinal  tensile  stress 

The  individual  fiber  breakage  occurs  when  the  applied  stresses  exceed  the  strength  of  the 
weakest  fibers  in  the  unidirectional  composite.  This  may  happen  as  the  service  loads  of  higher 
values  are  being  applied.  Initially,  such  fiber  failures  and  associated  intralaminar  voids  do  not 
significantly  affect  overall  strength  of  the  composite  as  the  fibers  failed  are  the  weakest  in  the 
plies  so  they  carry  a  relatively  small  amount  of  load.  Such  isolated  failures  are  usually 
noninteractive  when  they  occur. 

The  intralaminar  voids,  which  are  formed  after  a  few  individual  fibers  fail,  may  interact  with 
existing  or  newly  formed  microcracks  in  the  surrounding  material.  As  a  result,  these  microvoids 
may  grow  along  with  the  wing-shaped  macrocracks  on  both  sides.  The  growing  macrocracks  may 
reach  the  neighboring  fibers  and  introduce  interfacial  stress  concentration  that  may  lead  either  to 
the  fiber  failure  by  cleavage  or  to  the  interfacial  shear  failure  (Tanimoto  and  Amijima  1975; 
Talreja  1987;  Reifsnider  1991).  The  interfacial  debonding  promotes  the  macrocrack  bridging, 
formation  of  a  macrovoid,  and  the  tensile  failure  of  the  neighboring  fibers.  When  interaction  and 
the  consequent  coalescence  of  neighboring  macrovoids  begin,  the  final  failure  of  a  composite  is 
imminent. 

3.1.2  Damage  Caused  by  Off-Axis  Fatigue  Loading.  Under  off-axis  fatigue  loading,  the 
probability  of  fiber  breakage  will  rapidly  decrease  as  the  off-axis  orientation  angle  increases.  The 
predominant  damage  mechanisms  involve  matrix  cracking  and  interfacial  debonding  between  the 
matrix  and  reinforcing  fibers.  Any  crack  nucleated  at  the  fiber-matrix  interface  is  subjected  to  the 
normal  stress  in  transverse  direction  and  the  tangential  stress  along  the  fiber  direction.  These 
stresses  result  in  opening  and  sliding  crack  displacements,  respectively.  The  relative  magnitudes 
of  these  crack  displacements  depend  on  the  off-axis  angle  defining  the  orientation  of  fibers.  This 
also  applies  to  the  plies  having  ±45°  or  ±30°  orientation  angles. 

The  opening  mode  of  the  interfacial  crack  growth  under  fatigue  conditions  has  a  greater  effect 
on  the  degradation  of  macroscopic  mechanical  properties  of  unidirectional  composites.  This 
fracture  mode  increases  porosity  of  the  composite  and  leads  to  complete  local  debonding  of  fibers 
and  interfacial  microvoids.  The  voids  with  complete  circumferential  debonding  of  the  interface 
increase  effective  porosity  of  the  composite  and  lead  to  macrovoids.  Interaction  between 
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interfacial  voids  and  matrix  microcracking  may  result  in  the  wing-shaped  macrocracks  on  both 
sides  of  a  microvoid.  Such  voids  usually  become  macrovoids  via  the  damage  mechanisms 
described  in  the  previous  section. 

The  crack-opening  mode  affects  the  fatigue  limit  that  is  defined  by  the  limiting  applied  strain 
below  which  no  crack  growth  can  occur  (Talreja  1987).  The  fatigue  limit  will  be  the  lowest  under 
the  transverse  loading  normal  to  the  fiber  direction.  At  the  maximum  off-axis  angle,  the  opening 
mode  of  crack  growth  will  be  the  only  active  damage  mechanism.  The  mechanical  properties  of 
the  fiber-matrix  interface  or  the  matrix  material  properties  determine  the  minimum  strain  required 
for  the  initiation  of  transverse  fiber  debonding. 

3.2  Damage  in  Cross-Ply  Laminates.  The  fatigue  damage  mechanisms  in  cross-ply 
laminates  include  interfacial  debonding  of  transverse  fibers  (in  the  90°  plies),  transverse  matrix 
cracking,  individual  fiber  breakage,  interfacial  shear  failure,  and  interlaminar  fracture  (Owen  1974; 
Bailey  and  Parvizi  1981;  Berglund,  Varna,  and  Yuan  1992).  Broutman  and  Sahu  (1969)  have 
shown  that  debonding  of  the  transverse  fibers  constitutes  the  first  damage  mechanism  activated 
under  fatigue  loading  in  cross-ply  laminates.  The  resulting  transverse  interfacial  cracks  may 
interact  with  the  matrix  microcracking  near  the  fiber-matrix  interphase.  This  often  leads  to  the 
bridging  of  the  existing  transverse  cracks.  The  macrocracks  then  grow  toward  interlaminar 
interfaces  and  cause  stress  concentrations  in  the  interlaminar  region.  Interlaminar  delamination 
may  follow  the  crack  growth,  since  the  most  likely  initiation  sites  for  delamination  are  the 
transverse  cracks  in  the  matrix  (Crossman  et  al.  1980;  Korczynsky  and  Morley  1981;  Kim  and 
Mai  1991;  O'Brien  1993). 

The  interfacial  and  transverse  cracking  in  the  90°  plies  is  a  progressive  type  of  damage.  Their 
growth  has  probabilistic  nature  (Fukunaga  et  al.  1984).  The  developed  transverse  cracks 
approximately  span  the  thickness  and  width  direction  of  the  90°  layer.  The  transverse  cracking 
may  reach  a  saturation  state  (i.e.,  the  characteristic  damage  state)  before  the  final  failure  of  a 
composite.  Such  a  network  of  macrocracks  not  only  creates  stress  concentrations  and  the 
delamination  sites,  but  may  also  induce  the  fiber  breakage  in  adjacent  0°  plies.  Such  fracture 
events  have  progressive  character,  as  opposed  to  the  breakage  of  individual  weak  fibers  in  the  0° 
plies,  which  is  a  nonprogressive  type  of  failure.  The  weak-fiber  breakage  may  result  in  a  damage 
state  when  the  neighboring  microvoids  start  interacting  with  each  other,  matrix  microcracks  or  the 
transverse  macrocracks.  Such  interaction  is  characteristic  for  the  last  stage  in  the  fatigue  damage 
accumulation  before  the  final  failure. 
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It  should  be  noted  that  when  the  [0n/90m]s  laminates  are  being  loaded,  there  are  two  competing 
fracture  modes  during  the  initial  stage  of  loading.  The  aforementioned  matrix  and  interfacial 
cracking  in  transverse  direction  can  be  accompanied  by  the  free-edge  delamination.  It  is  caused 
by  the  interlaminar  tensile  stress  at  the  free  edge  (Pipes  and  Pagano  1970).  Wang  and  Crossman 
(1980)  demonstrated  that  the  free-edge  effects  are  usually  small  in  the  [0/90]s  type  composites 
and  that  edge  delamination  cannot  occur  alone.  Consequently,  Korczynsky  and  Morley  (1981), 
Crossman  and  Wang  (1982),  O'Brien  et  al.  (1982),  Kim  and  Mai  (1991),  and  O’Brien  (1993)  have 
been  studying  the  initiation  of  delamination  from  existing  transverse  cracks.  Transverse  cracking 
may  also  cause  local  longitudinal  splitting  (Harik  et  al.  1999). 

3.3  Damage  in  Woven  Fiber  Composites.  The  fatigue  damage  in  woven  fiber  composites 
begins  with  the  nucleation  of  numerous  cracks  during  the  first  cycle  of  loading  (Tanimoto  and 
Amijima  1975).  The  crack  density  then  gradually  reaches  a  constant  value  for  the  characteristic 
damage  state.  The  cracks  first  form  in  the  resin  area  near  the  fibers  perpendicular  to  the  load 
direction  (i.e.,  transverse  fibers).  Later,  the  cracks  grow  through  the  resin  matrix  surrounding  the 
transverse  fibers  toward  the  adjacent  longitudinal  fibers.  These  cracks  may  branch  off  when  they 
reach  any  material  inhomogeneity  such  as  a  fiber  or  fiber-matrix  interphase.  Depending  on  the 
interphasial  properties,  a  propagating  crack  may  be  “attracted”  or  “repelled”  by  the  encountered 
inhomogeneity  (Patton  and  Santare  1993). 

The  woven  laminates,  similar  to  the  commonly  used  glass/epoxy  and  graphite/epoxy  systems, 
usually  experience  matrix  cracking  under  tensile  stresses  greater  than  the  strength  of  plies.  During 
the  loading  process  of  these  relatively  brittle  matrices,  the  progressive  formation  of  cracks  results 
in  gradual  changes  in  the  compliance  of  the  composite.  The  progressive  nucleation  of  cracks 
continues  until  the  CDS  (i.e.,  characteristic  damage  state)  is  reached  in  the  crack-saturated  matrix 
(Reifsnider  1977).  The  final  crack  spacing  is  the  same  for  cyclic  and  fatigue  loading.  It  is  also 
independent  of  load  history,  environmental  conditions,  and  residual  and  moisture-induced 
stresses. 

It  is  important  to  know  whether  the  transverse  matrix  cracks  remain  in  the  resin  phase  or 
propagate  into  the  resin/interphase  areas  between  the  transverse  fibers.  If  cracks  remain  in  the 
resin  phase,  then  there  should  be  no  considerable  damage  and  corresponding  property  degradation 
in  the  composite  (Tanimoto  and  Amijima  1975).  The  matrix  cracks  may  also  form  parallel  to 
woven  fabric  lamina.  Their  length  and  density  increase  with  an  increasing  number  of  loading 
cycles.  These  cracks  may  also  reach  a  steady  state  as  well. 
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Damage  in  woven  composites  and  other  composite  systems  can  be  characterized  by  a  number 
of  methods.  The  commonly  used  S-N  curves  describe  the  residual  strength  of  composites,  which 
depends  on  the  number  of  cycles  to  failure.  Talreja  (1987)  proposed  the  strain-based  fatigue-life 
diagrams  that  characterize  the  residual  strain-to-failure  vs.  number  of  loading  cycles.  Talreja's 
fatigue-life  diagrams  provide  a  conceptually  useful  way  of  mapping  the  fatigue  damage 
mechanisms  onto  a  strain-vs.-cycle  diagram.  The  diagrams  do  not  characterize  any  of  the  rate¬ 
controlling  parameters  that  could  be  derived  only  from  more  quantitative  deterministic  or 
statistical  models  for  fatigue  damage  accumulation. 

4.  Modeling  of  Damage  Accumulation 

The  fatigue  damage  accumulation  models  should  reflect  the  physics  of  microscopic  damage 
initiation,  the  mechanics  of  crack  propagation  or  void  growth  within  composite  materials,  and  the 
changes  in  the  macroscopic  material  damage  state.  The  evolution  of  damage  in  laminated 
composites,  which  are  inhomogeneous  and  highly  anisotropic,  is  inherently  complicated.  The 
finite  element  methods  (Bathe  1982;  Zienkiewicz  and  Taylor  1989;  Kaliakin  1996)  proved  to  be 
powerful  tools  for  most  accurate  and  robust  modeling  of  nonlinear  deformation  and  fracture  of 
multiphase  materials  (Bogetti,  Gillespie,  and  Lamontia  1994;  Bogetti,  Hoppel,  and  Bums  1995; 
Harik  1997),  interfacial  fracture  (Needleman  1987;  Li  and  Kaliakin  1993;  Nath,  Fenner,  and 
Galiotis  1996;  Harik  1997),  evolution  of  voids  (Needleman  1972;  Becker  et  al.  1988)  and 
viscoelastic  effects  (Brinson  and  Knauss  1992;  Chen,  Davila,  and  Baker  1998).  However,  there 
are  many  difficulties  with  numerical  stability  and  error  analysis  of  nonlinear  finite  element 
simulations  (Nochetto  1990;  Szabo  and  Babuska  1991). 

The  patterns  of  fatigue  damage  accumulation  are  influenced  by  the  level  of  loads  relative  to 
the  ultimate  strength  of  composite  constituents,  the  rate  of  loading,  and  various  environmental 
conditions.  The  damage  accumulation  process  can  be  characterized  by  the  loss  of  stiffness, 
residual  strength,  residual  strain-to-failure  (Talreja  1987),  residual  fatigue  life  (Case  and 
Reifsnider  1998),  and  nondimensional  groups  of  physical  parameters  (Bridgman  1922).  The 
process  of  fatigue  damage  accumulation  can  be  divided  into  several  stages  such  as  initiation  of 
damage,  growth  of  pre-existing  defects,  noninteractive  and  interactive  evolution  of  multiple  cracks 
and  microvoids,  and  unstable  damage  accumulation  (Talreja  1987;  Reifsnider  1991;  Hahn  1979). 
Each  such  stage  can  be  described  by  one  or  more  material  damage  state  (Reifsnider  1994). 
Significant  loads  and  strain  rates  associated  with  the  LCF  may  have  considerable  effects  on  each 
stage  of  fatigue  damage  accumulation. 
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4.1  An  Initial  Damage  State  of  Composites.  The  initial  macroscopic  properties  of  PMCs 
depend  on  the  extent  of  processing-induced  damage  (e.g.,  microvoids,  cracking,  and  imperfect 
interfaces)  and  environmental  effects  that  constitute  an  initial  damage  state  (IDS)  of  a  composite 
(Harik  1997).  In  most  theoretical  analyses,  the  ideal  structural  state  (ISS)  is  assumed  as  the 
initial  state  of  a  composite  system  without  any  damage.  The  ISS  characteristics  include  a 
flawless  lay-up,  perfect  interfacial  and  interlaminar  bonding,  uniform  fiber  radii,  negligible  fiber 
waviness,  controlled  distributions  of  fibers,  absence  of  microvoids  and  flaws,  negligible  residual 
stresses,  and  uniform  degree  of  cure  (Hull  1981).  Variations  in  constituents  properties  and 
manufacturing  processes  may  result  in  different  IDSs  for  composites  having  the  same  ISS. 
Different  IDSs  result  in  the  structural  states  having  distinct  macroscopic  mechanical  properties. 

In  the  previous  section,  various  damage  evolution  patterns  are  discussed  for  several  types  of 
laminate  composites.  Most  of  these  patterns  involve  initiation  of  transverse  interfacial  cracks 
and/or  transverse  matrix  cracking.  Interfacial  defects  and  cracking  may  be  also  initiated  by 
thermal  loads  (Sottos  1990)  and  by  material  processing  (Harik  1997).  Such  cracking  defines  the 
initial  damage  state  of  a  composite.  The  stress  analyses,  which  are  based  on  the  unconstrained 
transverse  tensile  strength  of  the  90°  layer,  can  predict  the  onset  of  transverse  cracking  (Garret 
and  Bailey  1977;  Bailey,  Curtis,  and  Parvizi  1979;  and  Adams,  Bowles,  and  Herakovich  1986). 
However,  the  crack  initiation  strains  predicted  by  these  methods  are  unrealistically  small  in  the 
case  of  thin  transverse  layers.  In  practice,  the  constraining  effect  of  the  0°  layers  on  the 
transverse  crack  growth  in  the  90°  is  significant.  It  results  in  much  higher  strains  required  for  the 
onset  of  transverse  cracking. 

The  fracture  mechanics  approach,  which  is  based  on  an  energy  criterion,  is  able  to  describe 
the  thickness  effect  constraining  the  crack  growth.  This  method  postulates  that  a  microcrack  will 
form  when  the  released  energy  due  to  crack  propagation  is  greater  than  some  critical  value.  This 
value  is  called  the  critical  energy  release  rate.  This  method  is  proved  to  be  effective  for  predicting 
transverse  cracking  in  brittle  thermoset  composites.  In  the  case  of  tough  thermoplastic 
composites,  this  method  is  not  as  effective  (Berglund  et  al.  1992).  The  ductile  aspects  of  material 
behavior  under  finite  strain  rates  can  be  simulated  by  various  viscous  material  models  (Bingham 
1922;  Nadai  1950;  Harik  1997;  Barbat  et  al.  1997). 

Micromechanical  analysis  of  fatigue  damage  may  take  into  account  the  initial  variations  in  the 
microstructural  mechanical  properties.  Manders  et  al.  (1983)  and  Gao  and  Reifsnider  (1993) 
introduced  micromechanical  models  that  analyze  the  effect  of  statistical  distributions  for  the 
strength  of  fiber,  matrix,  and  the  fiber-matrix  interface.  The  activation  of  different  damage  and 
failure  modes  is  investigated  by  employing  a  number  of  criteria  for  the  onset  of  fracture  and 
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failure.  For  different  levels  of  loading,  the  degradation  of  mechanical  properties  and  stress 
redistributions  are  examined  for  both  monotonic  and  cyclic  loading.  These  models  can  be  used  to 
carry  out  parametric  studies  and  sensitivity  analyses  for  various  property  distributions  (Gao  and 
Reifsnider  1993).  As  a  result,  this  methodology  can  be  useful  for  structural  optimization  of 
composites. 

4.2  Micromechanical  Analysis  of  Individual  Cracks.  Microscopic  analysis  of  damage 
resolves  localized  stress  and  strain  concentrations  and  individual  fracture  events  (Love  1944; 
Timoshenko  and  Goodier  1951;  Landau  and  Lifshitz  1986).  Many  analytical  models  have  been 
proposed  to  solve  the  transverse  cracking  problem  (e.  g.,  Reifsnider  1977;  Parvizi  and  Bailey 
1978;  Bailey,  Curtis,  and  Parvizi  1979;  Flaggs  1985;  Dvorak  and  Laws  1987;  Gillespie  and 
Hansen  1996;  Akshantala  and  Talreja  1998).  Garret  and  Bailey  (1977)  used  the  shear-lag  theory 
to  derive  a  second-order  differential  equation  for  the  amount  of  stress,  Ag,  transferred  from  the 
0°  ply  to  the  90°  damaged  ply: 


d1  Ag 

~w 


+  02Ag  =  0, 


(1) 


where  b,  —  xh  and  the  constant  0  is  given  by  02  =  G^Q  where  Cx  =  —  +  — — .  Here,  Gjj,  ET, 

ET  AEt 

and  EL  are  the  shear,  transverse,  and  longitudinal  moduli,  respectively,  for  the  unidirectional 
composite.  In  the  derivation  of  the  shear-stress  transfer  coefficient,  0,  the  through-thickness 
uniformity  of  the  displacement  in  the  transverse  plies  in  x-direction  was  assumed. 


A  number  of  researchers  made  various  improvements  to  the  Garret-Bailey  analysis.  Manders 
et  al.  (1983)  extended  the  Garret-Bailey  analysis  to  include  the  effects  of  neighboring 
microcracks.  Ogin  and  Smith  (1987)  assumed  that  the  displacement  in  the  90°  plies  is  parabolic 
in  z,  and  derived  the  Garret-Bailey  equation  with 


02  =  3GttCv  (2) 

Han,  Hahn,  and  Croman  (1987)  obtained  die  same  results.  Reifsnider  (1977)  and  Dvorak  and 
Laws  (1987)  introduced  a  shear  transfer  layer  between  the  ply  groups,  which  is  characterized  by 
the  effective  shear  stiffness,  0,  as  an  adjustable  parameter.  Reifsnider's  analysis  yields 
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where  G  is  the  shear  modulus  of  the  shear  transfer  layer  and  to  is  its  thickness.  In  this  approach, 
the  shear  stiffness,  G/to,  of  the  shear  stress  transfer  layer  is  an  unknown  parameter  that  must  be 
determined  by  fitting  experimental  data. 

The  shear-lag  method  is  simple  and  gives  reasonable  predictions  of  stiffness  reductions. 
However,  it  ignores  the  shear  and  through-thickness  deformation  due  to  the  opening  displacement 
of  a  transverse  crack.  Moreover,  the  shear-lag  solution  predicts  a  non-zero  shear  stress  on  the 
transverse  crack  surfaces,  which  violates  the  no-shear  boundary  condition.  Since  the  shear-lag 
approach  is  essentially  one-dimensional,  it  can  not  provide  complete  stress  distributions. 

Hashin  (1985,  1986)  proposed  a  two-dimensional  analysis  of  the  mechanical  response  in  the 
x-z  plane  by  using  the  principle  of  minimum  complementary  energy  in  the  variational  framework. 
In  contrast  to  the  shear-lag  models,  Hashin's  closed-form  solution  satisfies  the  no-shear  stress 
boundary  condition.  It  also  estimates  the  interlaminar  stresses.  The  through-thickness  variations 
of  the  axial  normal  stresses  are  assumed  to  be  small.  Varna  and  Berglund  (1991,  1992)  modified 
Hashin's  variational  model  to  take  into  account  a  normal-stress  gradient  in  the  0°  layer.  However, 
the  stress  gradient  in  the  90°  layer  was  still  neglected.  All  of  the  aforementioned  models  describe 
stress  distributions  around  a  transverse  crack  that  spans  the  whole  thickness  and  the  whole  width 
of  the  transverse  layer.  The  growth  process  of  the  transverse  crack  was  not  considered.  Wang 
and  Crossman  (1980)  suggested  an  energy  method  for  the  investigation  of  initiation  and  growth  of 
transverse  cracks  and  edge  delamination  in  composite  laminates.  They  assumed  that  a  microcrack 
exists  in  the  90°  layer  of  the  matrix.  Conditions  for  the  stable  crack  growth  were  investigated  in 
the  framework  of  "effective  flaws." 

Akshantala  and  Talreja  (1998)  proposed  a  mechanistic  model  of  the  evolution  of  transverse 
cracking  in  cross-ply  laminates  that  are  subjected  to  cyclic  tension  in  the  longitudinal  direction. 
The  unique  feature  of  their  model  is  that  it  takes  into  account  delamination  associated  with 
transverse  cracks  so  the  progressive  delamination  induces  further  formation  of  transverse  cracks. 
In  the  region  between  transverse  cracks,  the  stresses  are  estimated  by  a  variational  approach, 
which  was  shown  to  yield  an  accurate  solution  away  from  the  crack  planes.  This  model  allows  an 
effective  prediction  of  the  transverse  crack  density  and  changes  in  the  macroscopic  elastic  moduli. 

4.3  Analysis  of  Delamination  Growth  in  PMCs.  Composite  laminates  develop  significant 
interlaminar  stresses  under  axial  tension  (Pipes  and  Pagano  1991).  Wang  and  Crossman  (1980) 
pointed  out  that  the  free-edge  effects  are  usually  small  in  the  composites  having  [0/90]s  lay-up, 
so  the  edge  delamination  cannot  occur  alone.  Hence,  the  combined  effect  of  matrix  cracking  and 
free  edge  delamination  should  be  investigated  (Salpekar  and  O'Brien  1991).  O'Brien  et  al.  (1982) 
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have  developed  an  effective  test  for  experimental  characterization  of  the  interlaminar  fracture 
toughness  of  composites  and  studies  of  the  interlaminar  crack  growth.  Whitcomb  (1991)  and 
O'Brien  (1993)  showed  that  the  initiation  of  delamination  occurs  at  existing  transverse  cracks  or 
other  sources  of  material  nonhomogeneity.  Inclusions  may  attract  or  repel  a  propagating  crack  so 
that  it  may  change  its  path  (Patton  and  Santare  1993). 

In  the  fatigue  performance  analysis  (e.g.,  MRLifel  1,  see  Appendix),  delamination  can  be  also 
assumed  to  be  caused  by  the  interlaminar  stresses  at  a  free  edge.  The  strain  energy  release  rate, 
G,  for  the  interlaminar  crack  is  then  estimated  by  a  fracture  mechanics  approach  in  conjunction 
with  the  laminated  plate  theory  (O'Brien  1991): 

G  =  ^Ekm-E-),  (4) 

where  £  is  the  strain,  t  is  the  thickness  of  the  laminate,  and  E[am  and  £*  are  laminate  moduli  before 
delamination  and  after  total  delamination,  respectively.  The  estimate  is  independent  of  the 
delamination  length,  but  it  takes  into  account  the  thickness  of  the  laminate.  This  results  in  the 
O'Brien-Paris  law. 


-j-  =  AGb  ,  (5) 

dn 

for  the  interlaminar  crack  growth  from  a  delaminated  edge  under  fatigue  loading.  Here,  a  is  the 
delamination  length,  n  is  the  number  of  cycles,  and  A  and  B  are  experimentally  determined 
constants.  This  model  for  the  interlaminar  fatigue  crack  growth  is  similar  to  the  Paris  law  for 
fatigue  crack  growth  in  metals. 

4.4  Micromechanical  Analysis  of  Voids.  Interfacial  debonding,  sliding  and  separation,  and 
interfacial  void  formation  are  among  interfacial  problems  that  may  occur  during  processing  and 
finite-strain-rate  LCF  loading.  These  interfacial  problems  can  later  lead  to  premature  failure. 
Recently  these  interfacial  problems  have  been  studied  by  a  number  of  researchers  (e.g.,  Xia  et  al. 
1994;  Budiansky,  Evans,  and  Hutchinson  1995;  Jasiuk  and  Kouider  1993;  and  Harik  and 
Caimcross  1999).  When  such  materials  are  subjected  to  compressive  or  extensional  loads, 
interfacial  voids  may  occur  by  further  decohesion  of  the  matrix  material  from  the  inclusions 
(Hashin  1991). 

Budiansky,  Hutchinson,  and  Slutsky  (1982)  investigated  the  evolution  of  isolated  spherical 
voids  in  an  infinite  linear  viscous  solid  subjected  to  various  bi-axial  stresses.  The  effect  of 
outside  pressure  on  the  final  shape  of  such  voids  was  extensively  analyzed.  Deformation  of 
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spherical  cavities  has  been  also  studied  by  Rice,  Rudnicki,  and  Simons  (1978)  for  a  class  of  fluid- 
infiltrated  elastic  materials.  For  viscous  materials,  Budiansky,  Hutchinson,  and  Slutsky  (1982) 
showed  that  tension  or  transverse  compression  loads  lead  to  elongated  ellipsoidal  voids.  Harik 
and  Caimcross  (1999)  showed  that  similar  tendencies  develop  in  plane  compression  flows  around 
cylindrical  inclusions,  although  the  deformation  of  cavities  is  no  longer  homogeneous.  Needleman 
(1987)  and  Lee  and  Batt  (1989)  studied  formation  of  interfacial  separation  and  evolution  of 
interfacial  voids  at  rigid  inclusions  in  an  elastic- viscoplastic  matrix.  In  particular,  it  was  noted 
that  the  shear  stiffness  parameter  of  the  phenomenologically  described  interface  had  insignificant 
effect  on  the  voids  studied.  In  the  case  of  compression  flows  around  cylindrical  particles  with 
perfectly  weak  interfaces,  Harik  and  Caimcross  (1999)  showed  that  the  interfacial  sliding  has 
rather  small  influence  on  the  voids  as  they  are  formed  by  predominantly  normal  interfacial 
separation.  Such  voids  may  grow  and  collapse  (Lee  and  Dawson  1993;  Lee  and  Mear  1994).  The 
monotonically  growing  voids  may  lead  to  void  coalescence  (Koplik  and  Needleman  1988)  and 
unstable  failure. 

4.5  Continuum  Damage  Mechanics  of  Multiple  Cracks  and  Voids.  In  contrast  to  the 
traditional  (differential)  fracture  mechanics,  the  Continuum  Damage  Mechanics  (CDM)  describes 
the  accumulation  of  multiple  cracking  or  voidage  and  their  effect  on  the  degradation  of 
macroscopic  mechanical  properties  (Kachanov  1958, 1986;  Chaboche  1981;  Lemaitre  1984,  1992; 
Krajcinovic  1989;  Talreja  1990;  Voyiadjis  1998,  1999).  Kachanov  (1958)  was  the  first  to 
describe  the  effect  of  multiple  voids  on  high-temperature  creep  behavior  by  introducing  an 
internal  damage  variable  D.  Macroscopic  effects  of  damage  are  described  by 


Y 
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and 


D  = 
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o0(l-D) 


(7) 


where  n,  p,  e0,  and  <r0  are  experimentally  determined  parameters  and  a  and  e  denote 
unidirectional  stress  and  strain,  respectively.  The  damage  variable  D  is  assumed  to  be  a  scalar 
function.  Although  there  is  contrary  microscopic  evidence  indicating  the  directional 
characteristics  of  damage,  the  scalar  damage  is  a  useful  modeling  variable  (Kachanov  1958,  1986). 
This  approach  is  relatively  simple,  and  it  can  be  readily  related  to  the  experimental  data  on 
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macroscopic  mechanical  properties.  More  general  models  were  developed  by  Lemaitre  (1984, 
1992),  Talreja  (1987,  1990,  1996),  Krajcinovic  (1989),  Arnold  and  Kruch  (1991),  Voyiadjis  and 
Echle  (1998),  and  Voyiadjis  and  Park  (1999). 

The  damage  variable  D  and  the  corresponding  accumulation  of  damage  can  be  characterized  by 
the  stiffness  loss  measured  during  testing  (Burr,  Hild,  and  Leckie  1995).  In  1984,  Ashby  and 
Dyson*  established  a  connection  between  the  damage  parameter  D  and  the  property  deterioration 
occurring  in  metals  subjected  to  a  high-temperature  environment.  A  functional  map  linking  the 
operating  conditions  and  various  damage  mechanisms  was  constructed.  In  1987,  Cocks  and 
Leckie*  recast  these  results  into  the  CDM  format.  In  1991,  Hall  and  Hayhurst*  demonstrated 
how  the  CDM  method  can  accurately  predict  the  type  and  growth  of  internal  damage  in 
engineering  structures  under  high-temperature  conditions.  Arnold  and  Wilt  (1993),  Voyiadjis  and 
Echle  (1998)  developed  computational  algorithms  based  on  CDM  models  for  deformation, 
damage  accumulation,  and  lifetime  prediction  of  composites. 

In  the  CDM  models,  the  matrix  constitutive  properties  are  locally  averaged  and  determined 
with  respect  to  the  damaged  local  volume.  The  effect  of  transverse  cracks  is  reflected  in  the 
constitutive  equation  via  introduction  of  homogenized  effective  layers  instead  of  the  cracked 
transverse  layers.  These  models  usually  provide  only  in-plane  stress  components  because  the 
stresses  are  derived  from  the  effective  constitutive  equations  and  classical  laminate  theory.  As  a 
result,  the  out-of-plane  stress  and  interlaminar  stresses  cannot  be  determined  by  the  CDM 
methods.  This  restricts  the  effectiveness  of  the  CDM  methods,  since  the  interlaminar  fracture 
constitutes  a  very  common  failure  mechanism  that  has  to  be  taken  into  account. 


5.  Fatigue  Life  Prediction 

Advanced  composite  materials  perform  well  in  weight-critical  structural  applications.  To 
ensure  the  weight  advantage  of  laminates,  one  needs  to  improve  their  fatigue  performance,  as  they 
are  prone  to  delamination  failures  under  fatigue  conditions  (O'Brien  1991).  This  is  especially 
important  for  Army  land  combat  systems  that  are  often  subjected  to  high  loads  and  moderate 
strain  rates  under  the  LCF  conditions.  In  order  to  improve  the  mechanical  behavior  of  these 
composites  under  any  fatigue  loading,  it  is  necessary  to  link  their  material  and  structural 
characteristics  to  their  fatigue  life  span.  Schaff  and  Davidson  (1997a,  b)  point  out  that  current 


*  See  Burr,  Hild,  and  Leckie  (1995). 
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life  prediction  models  are  ineffective  for  the  PMCs.  Counterbalancing  of  this  deficiency  leads  to 
“overdesigned”  structures  and  large  factors  of  safety,  which  inevitably  result  in  heavier  and  more 
costly  structural  components.  Hence,  heavier  vehicles  have  much  lower  vehicle  performance, 
which  translates  into  reduced  mobility. 

The  development  of  service  life  prediction  models  requires  not  only  understanding  of  the 
basic  mechanisms  of  material  degradation  but  also  reliable  models  of  damage  accumulation  and 
damage  effects  on  the  macroscopic  mechanical  properties  (Sendeckyj  1990).  The  presence  of 
high  loads  and  significant  strain  rates  under  LCF  loading  conditions  requires  nonlinear  modeling 
capabilities.  The  fatigue  life  predictions  can  be  based  on  either  the  rate  of  damage  accumulation 
or  on  the  rate  of  property  degradation.  These  approaches  can  be  characterized  as 
micromechanical  or  mechanistic  and  phenomenological  or  macromechanical,  respectively. 
Micromechanical  models  quantitatively  account  for  the  microstructural  effects  and  progression  of 
damage.  Macroscopic  phenomenological  models  are  based  on  macroscopic  properties  (strength, 
stiffness,  etc).  A  hybrid  approach  involving  prediction  of  damage  growth  and  property 
degradation  can  also  be  used. 

5.1  Micromechanical  Fatigue  Models.  Micromechanical  fatigue  models  provide  a 
quantitative  account  for  the  effects  of  a  particular  microstructure  and  the  transient  evolution  of 
microscopic  damage  in  composites.  Such  mechanistic  models  are  independent  of  lay-ups  and 
type  of  loading,  so  they  are  adaptable  to  various  geometric  variations.  These  models  usually 
require  minimal  experimentally  obtained  input.  As  a  result,  such  models  promise  to  be  useful  for 
wide  variety  of  existing  and  new  composite  materials.  However,  the  complex  nature  of  fatigue 
phenomena  in  highly  anisotropic  materials  poses  major  challenges  for  understanding  and  accurate 
modeling  of  the  physical  processes  occurring  on  the  microscopic  level. 

Gao  and  Reifsnider  (1993)  introduced  micromechanical  models  that  analyze  the  effect  of 
statistical  distributions  for  the  strength  of  fiber,  matrix,  and  the  fiber-matrix  interphase  on  the 
macroscopic  mechanical  properties.  The  activation  of  different  failure  modes  is  investigated  by 
employing  a  number  of  failure  criteria.  For  different  levels  of  loading,  the  degradation  of 
mechanical  properties  and  stress  redistributions  are  examined  for  both  static  and  cyclic  loading. 
Dzenis,  Joshi,  and  Bogdanovich  (1993)  and  Dzenis  (1996)  proposed  other  micromechanical 
damage  models  which  may  include  stochastic  effects  in  damage  evolution.  These  models  can  be 
used  to  carry  out  parametric  studies  and  sensitivity  analyses  for  various  property  distributions 
(Gao  and  Reifsnider  1992).  As  a  result,  this  methodology  can  be  useful  for  structural 
optimization  of  composites. 
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Talreja  (1990)  suggested  a  tensorial  representation  of  various  damages  within  a  composite 
material.  A  given  damage  entity  within  volume,  V,  is  assumed  to  be  bounded  by  a  surface,  S, 
with  a  uniquely  defined  unit  vector,  n.  The  mechanical  influence  of  each  point  on  the  volume,  V, 
is  described  by  a  vector,  a.  As  a  result,  a  second  order  tensor,  d,  can  be  defined  by 

dy  J  apjdS  (8) 

s 

If  N  damage  modes  are  present,  then  for  each  a  th  mode,  one  may  define  the  damage  tensor 


where  a=  1,2,  ...  N  and  ka  is  the  number  of  damage  entities  in  the  ath  mode.  The  damage 
influence  vector,  a,  can  be  decomposed  into  its  normal  and  tangential  components.  The  normal 
component  is  of  special  interest  for  the  transverse  matrix  cracking.  For  small  strains  and  low 
concentrations  of  damage,  the  stiffness  of  the  composite  material  can  be  easily  related  to  the 
damage  tensor  Dy  (Talreja  1990).  Therefore,  the  presence  of  damage  affects  the  initial  symmetry 
properties  of  the  material  and  its  stiffness.  Therefore,  the  degradation  of  material  stiffness  can 
be  predicted. 

5.2  Phenomenological  Fatigue  Models.  Phenomenological  models  are  effective  simulation 
tools  for  analyzing  macroscopic  behavior  of  composites.  The  structure  of  these  models  may 
vary  from  simple  empirical  or  semi-empirical  rules  to  complex  concepts  rooted  in  the  continuum 
damage  mechanics  (Reifsnider  1991;  Talreja  1994;  Case  and  Reifsnider  1998;  Chamis  1999;  Tsai 
1999).  Macroscopic  models  are  especially  promising  for  specific  industrial  applications  where 
achieving  short-term  modeling  goals  is  important.  The  simpler  models  tend  to  be  reliable  only  for 
a  narrow  group  of  composite  materials  that  have  many  structural  similarities.  This  limitation 
stems  from  the  fact  that  the  phenomenological  models  do  not  quantitatively  account  for  evolution 
of  damage  in  composites.  The  more  differences  between  the  material’s  microstructure  exist,  the 
more  variety  in  the  patterns  of  damage  evolution.  Therefore,  these  structural  differences  result  in 
the  variations  in  the  damage  accumulation  effects  on  the  macroscopic  behavior  of  composites. 
The  macroscopic  effects  of  typical  microstructures  can  be  characterized  by  comprehensive 
fatigue  testing. 

The  structure  of  composite  materials  may  also  vary  on  the  ply  level.  Therefore,  the  stacking 
sequence  and  other  characteristics  of  lay-ups  may  affect  the  macroscopic  properties  of 
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composites  as  well.  The  microstructural  effects  also  include  local  stress  distributions  and  stress 
concentrations  that  vary  from  one  composite  to  another,  even  if  they  have  the  same 
microstructure.  Such  variations  result  from  pre-existing  or  processing-induced  defects,  which  are 
also  not  accounted  for  in  phenomenological  models.  Local  stress  distributions  may  also  change  as 
the  characteristics  of  loading  change.  Hence,  simple  macroscopic  models  are  sensitive  to  the  type 
of  fatigue  loading.  This  sensitivity  of  empirical  models  to  the  initial  microstructure  is  typical  for 
mathematical  modeling  of  any  ill-conditioned  physical  problem.  The  effects  of  typical  structural 
variations  can  be  often  evaluated  by  additional  experimental  testing.  Multiple  experiments  are 
also  required  for  fatigue  analyses  based  on  the  time-temperature  superposition  hypothesis  (Tsai 
1999).  This  dependency  on  large  amounts  of  experimental  input  for  each  type  of  material,  lay¬ 
up,  and  loading  is  a  major  disadvantage  of  all  phenomenological  models. 

5.2.1  Models  Based  on  Strength  and  Stiffness.  The  phenomenological  life  prediction  models 
characterize  the  degradation  of  macroscopic  mechanical  properties,  such  as  strength,  stiffness, 
etc.  The  material's  strength  and  stiffness  are  the  primary  characteristics  of  mechanical  behavior 
that  could  be  easily  monitored.  As  a  consequence,  there  are  stiffness-based  and  strength-based 
models  for  fatigue  life  predictions.  The  fatigue  failure  of  a  composite  occurs  when  the  current 
stress  applied  is  equal  to  (or  greater  than)  its  residual  strength.  It  is  physically  natural  to  employ 
this  failure  criterion  in  the  strength-based  fatigue  models  (Case  and  Reifsnider  1998;  Appendix). 
These  models  are  often  characterized  as  “wearout”  models. 

The  models,  which  use  the  composite's  stiffness  as  their  primary  variable,  additionally 
require  formulation  of  fundamental  failure  criteria.  Hahn  and  Kim,  O'Brien  and  Reifsnider, 
Whitworth  and  Farrow  (Schaff  and  Davidson  1997a)  introduced  failure  criteria  based  on  the 
secant  modulus.  The  failure  criteria  based  on  the  static  strain  to  failure  were  proposed  by  Hwang 
and  Han  (1986),  and  Poursartip,  Ashby  and  Beaumont  (1986).  The  stiffness-based  fatigue 
models  associated  with  the  aforementioned  criteria  provide  reasonable  fatigue  life  predictions  for 
constant  amplitude  and/or  two-stress  amplitude  loadings.  The  degradation  rates  for  the  Young's 
modulus  and  the  strength  of  composites  can  be  evaluated  via  experimental  testing.  The 
experimentally  determined  rates  can  be  characterized  by  such  concepts  as  the  fatigue  degradation 
modulus  (Hwang  and  Han  1986)  and  a  factor  for  degradation  of  strength  (Schaff  and  Davidson 
1997a,  b).  These  fatigue  parameters  are  often  used  in  the  phenomenological  life  prediction 
models  to  characterize  the  degradation  of  mechanical  properties  or  the  residual  strength,  the 
remaining  fatigue  life,  etc.  It  should  be  emphasized  that  the  rates  of  material  property 
degradation  are  important  material  parameters  that  should  be  known  during  any  design  process. 
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For  instance,  initially  stiff  composites  may  have  worse  fatigue  performance  than  some  more 
compliant  composite  materials  because  of  their  high  rates  of  the  modulus  degradation. 

5.2.2  Statistical  Fatigue  Models.  The  residual  strength  and  fatigue  life  of  composites  are 
statistical  quantities,  in  general  (Fukunaga  et  al.  1984;  Talreja  1987;  Chamis  1999).  Two- 
parameter  Weibull  functions  are  commonly  used  in  the  strength-based  fatigue  models  to  describe 
the  residual  strength  distribution  after  arbitrary  load  history.  These  functions  can  also  describe 
the  probability  of  the  composite's  failure  and  fatigue  life  distributions  after  an  arbitrary  number 
of  fatigue  cycles.  The  two-parameter  Weibull  functions  are  defined  by  the  scale,  which 
represents  the  63.2  percentile  of  the  distribution,  and  the  shape,  which  characterizes  the  degree  of 
scatter  in  the  statistical  data.  Both  the  scale  and  shape  parameters  can  be  determined  by 
experimental  testing  and  the  method  of  maximum  likelihood  (Talreja  1987). 

The  shape  parameter  for  strength  at  zero  cycles,  5/0),  must  equal  the  static  shape  factor,  Bs. 
The  shape  parameter  usually  decreases  with  increasing  cycling,  as  the  residual  strength 
distributions  become  more  disperse  during  fatigue  testing  (Schaff  and  Davidson  1997a).  The 
range  of  values  involved  in  the  residual  strength  distribution  becomes  wider  as  the  number  of 
fatigue  cycles  increases.  The  Weibull  scale  parameters  are  defined  as  the  63.2  percentile  of  the 
respective  distribution  functions.  Experimental  data  on  the  initial  static  strength  and  the  residual 
strength  of  composites  can  be  represented  by  the  Weibull  scale  parameters.  Schaff  and  Davidson 
(1997a,  b)  developed  models  for  constant  amplitude  and  two-stress  amplitude  fatigue  with  a 
reduced  number  of  experimentally  determined  parameters.  The  Schaff-Davidson  models  are 
based  on  the  following  assumptions: 

•  environmental  and  frequency-related  effects  are  negligible, 

•  the  residual  strength,  R(ri),  initially  equals  the  static  strength,  Ro  (it  is  also  assumed  to  be  a 
monotonically  decreasing  function  of  fatigue  cycles,  n,n>N), 

•  the  residual  strength,  R(n),  the  static  strength,  R0,  and  the  final  cycle  number,  N,  are 
assumed  to  be  the  Weibull  scale  parameters  (i.e.,  they  are  defined  as  the  63.2  percentile  of 
their  respective  distribution  functions), 

•  the  Weibull  shape  parameter,  Bf,  is  a  linear  decreasing  function  of  fatigue  cycles  (i.e.,  the 
residual  strength  distribution  range  of  values  becomes  monotonically  wider  during 
fatigue), 

•  the  failure  occurs  when  the  residual  strength,  R(n),  equals  the  peak  stress,  Smax, 
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•  the  stress  ratio,  rs  =  is  constant  during  fatigue,  and 

•  the  residual  strength  relation  is  given  by  a  power  law,  which  is  defined  by  the  rate  of 
strength  1oss,/(Rq,  rs,  Smax),  and  the  strength  degradation  parameter,  v. 

The  strength  degradation  parameter  v  can  characterize  a  wide  variety  of  the  degradation  rates 
of  material  properties  under  fatigue  loading.  In  the  case  when  v  <  1,  there  is  a  rapid  degradation 
of  strength  in  the  beginning  of  service  life.  Linear  strength  degradation  corresponds  to  v  =  l .  The 
case  when  v  »1  is  characterized  by  the  “sudden  death”  behavior.  The  fatigue  life  predictions  by 
the  Schaff-Davidson  model  compare  well  with  experimental  results  for  low-high  and  high-low 
two-stress  amplitude  fatigue  tests  (Schaff  and  Davidson  1997b).  The  “cycle  mix”  effect  of  the 
changing  loading  sequences  has  been  taken  into  account  by  introducing  a  “cycle  mix  factor” 
(Schaff  and  Davidson  1997b).  It  is  a  scale  parameter  used  for  the  degradation  of  the  residual 
strength  during  the  transition  cycles.  As  a  result,  the  model  can  be  also  used  to  simulate  the  pulse 
loading  conditions.  Probabilistic  sensitivity  factors  can  be  also  used  to  account  for  uncertainty  in 
the  performance  and  durability  evaluation  of  composite  structures  (Chamis  1999). 


6.  Recommendations  for  ARL  LCF  Program 

The  LCF  conditions  are  unique  to  many  Army  land  combat  systems  and  are  not  well 
understood.  Army  engineers  need  to  develop  a  clear  physical  understanding  of  LCF  and  the 
effects  of  material  microstructure  on  fatigue  damage  processes,  a  methodology  for  LCF 
characterization  of  PMCs,  and  novel  LCF  models  for  optimization  of  designs  and  realistic  service 
life  prediction  of  Army  combat  systems  (e.g.,  gun  components,  integral  armor,  rotorcraft 
applications,  etc.).  The  broad  range  of  fatigue  problems  encountered  in  Army  structural 
applications  (see  section  2)  can  be  addressed  only  by  a  comprehensive  research  and  development 
program.  Significant  organizational  and  research  efforts  should  be  focused  on  establishing  a 
knowledge  base  and  infrastructure  for  a  coordinated  ARL  LCF  Program.  Such  a  program  will 
enable  ARL  engineers  to  meet  current  and  future  Army  needs  for  solutions  to  numerous  fatigue 
related  problems  in  design  and  repair  of  composite  structures. 

The  conceptual  framework  for  the  ARL  LCF  Program  has  to  reflect  the  Army  needs  for  the 
design-for-fatigue  methodology  and  robust  predictive  capabilities  for  microstructural  design 
optimization  and  LCF  life  prediction  of  structures  which  have  thick  sections  and  complex 
anisotropic  microstructure.  In  order  to  develop  reliable  LCF  models,  the  effects  of  fiber-matrix 
interphase  and  adhesive  joints  in  hybrid  composites  have  to  be  accounted  for.  These  issues  are 
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unique  to  mechanics  of  composite  materials  and  have  not  been  examined  in  the  studies  dealing  with 
fatigue  of  metals.  The  effects  of  chemical  and  hydro-thermal  environment  on  material  degradation 
processes  and  accumulation  of  damage  also  have  unique  features  which  are  specific  to  composite 
materials  used  by  the  Army. 

6.1  Recommendations  for  Design-for-Fatigue  Research.  Section  2.1  includes  some 
recommendations  for  future  research  concerning  design  methodologies,  which  can  be  summarized 
as  follows. 

•  Assess  the  structural  stress  analysis  data  on  specific  Army  composite  structures  in  order 
to  identify  the  critical  structural  elements,  the  types  of  critical  loads,  and  expected  damage 
and  failure  modes. 

•  Carry  out  experimental  fatigue  testing  of  specific  composite  materials  used  in  Army  land 
combat  systems.  Such  experimental  characterization  would  yield  S-N  curves,  residual 
strength  data,  stiffness  degradation  rates,  etc. 

•  Develop  methodologies  for  fatigue  life  assessment  based  on  either  the  experimentally 
determined  property  degradation  rates  or  phenomenological  fatigue  life  prediction  models 
involving  initial  damage  criteria  or  one-cycle  fracture  analysis. 

•  Implement  phenomenological  fatigue  models  into  the  LAMPAT's  structural  analysis 
capabilities.  Develop  design  methodologies  that  would  include  an  assessment  of  fatigue 
performance  based  on  engineering  fatigue  analysis. 

•  Develop  new  phenomenological  and  micromechanical  fatigue  modeling  capabilities  that 
can  account  for  microstructural  effects  such  as  interphasial  effects,  interfacial  damage, 
stacking  sequence  effects,  etc. 

•  Evaluate  current  state-of-the-art  fatigue  simulation  technology  such  as  MRLife  simulation 
code  (see  Appendix),  which  is  suited  for  fatigue  performance  analysis  of  a  wide  variety  of 
polymer  and  ceramic  composite  systems.  The  use  of  MRLife  code  would  accelerate  the 
progress  in  the  design-for-fatigue  efforts. 

•  Assess  the  possibility  of  linking  the  modeling  capabilities  of  MRLife  fatigue  performance 
simulation  code  with  the  structural  analysis  capabilities  of  commercial  FEA  codes  (e.  g., 
ANSYS,  ABAQUS,  etc.)  and  with  composites-specific  postprocessing  software 
programs  like  LAMP  AT. 
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6.2  Recommendations  for  the  LCF  Characterization.  Section  2.2  includes  some 
recommendations  for  future  research  concerning  LCF  for  the  ARL  LCF  Program,  which  can  be 
summarized  as  follows. 

•  Assess  the  effect  of  high  loads  and  various  strain  rates  on  the  transient  stress 
distributions  in  critical  structural  elements  of  Army  composite  structures.  Characterize 
the  associated  nonlinear  material  behavior  and  identify  expected  damage  and  failure 
modes. 

•  Perform  experimental  fatigue  testing  of  test  coupons  representing  the  critical  structural 
elements  and  carry  out  detailed  nondestructive  monitoring  of  fatigue  damage  accumulation 
processes  for  cyclic  loading  above  50%  of  the  ultimate  strength. 

•  Develop  conceptual  maps  for  physical  understanding  of  the  mechanisms  and  phenomena 
associated  with  fatigue-induced  failure  in  PMCs  under  LCF  conditions. 

•  Evaluate  the  effects  of  LCF  conditions  on  typical  damage  mechanisms  and  the  kinetics  of 
damage  accumulation  processes. 

•  Examine  evolution  of  ballistic  damage  under  LCF  conditions,  develop  models  for  LCF 
assessment  of  ballistically  damaged  PMC  structures,  and  validate  these  models. 

•  Develop  novel  LCF/PVF  models  accounting  for  fiber-matrix  interphase  effects  and 
realistic  service  life  prediction  methodology  for  various  Army  composite  systems. 

•  Develop  a  design  optimization  methodology  based  on  service  life  prediction  capabilities 
for  Army  structures  under  LCF  conditions. 

6.3  Recommendations  for  Investigation  of  Thick-Section  Composites.  The  thick- 
section  PMCs  are  unique  to  Army  systems  and  should  serve  as  a  focal  point  for  the  ARL  LCF 
Program.  Section  2.3  includes  some  recommendations  for  future  studies  of  thick-section 
composites,  which  can  be  summarized  as  follows. 

•  Identify  critical  thick-section  structural  elements  in  Army  land  combat  systems  and 
assess  the  types  of  critical  loads  and  expected  damage  and  failure  modes.  This 
assessment  should  be  based  on  the  structural  stress  analysis  data  for  specific  Army 
composite  structures. 
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•  Identify  similarities  and  differences  in  mechanical  behavior  of  the  thick-section  structures 
in  Army  applications  and  other  thick-section  homogeneous  and  heterogeneous  structures 
used  in  industry.  This  analysis  should  be  based  on  existing  theoretical  and  experimental 
data. 

•  Select  appropriate  test  fixtures  that  could  be  scaled  up  and  used  for  comprehensive 
fatigue  characterization  of  the  test  coupons  representing  the  critical  thick-section 
structural  elements  of  interest  to  Army. 

•  Identify  nondestructive  damage  evaluation  techniques  (e.g.,  x-rays,  computed 
tomography,  acoustic  wave  scattering,  infrared  thermal  imaging,  and  embedded  fiber 
optics),  which  would  be  appropriate  for  detailed  nondestructive  monitoring  of  fatigue 
damage  accumulation  in  the  thick-section  composites. 

•  Characterize  different  types  of  fatigue  damage,  and  their  unique  features,  and  determine 
the  dominant  fatigue  damage  mechanisms  in  the  thick-section  composites. 

•  Develop  various  "damage  healing"  techniques  (e.g.,  heating  of  magnetic  microscopic 
particles  under  loading)  for  different  types  of  fatigue  damage  in  the  thick-section 
composites,  especially  for  interlaminar  delamination. 

•  Develop  conceptual  maps  describing  the  accumulation  processes  of  fatigue  damage  in  the 
thick-section  composites  under  a  wide  range  of  loading  conditions. 

•  Identify  phenomenological  fatigue  models  that  can  predict  the  stiffness  degradation  rates 
and  fatigue  life  of  the  thick-section  composites. 

•  Develop  design  optimization  methodologies  that  would  include  an  assessment  of  fatigue 
performance  based  on  engineering  fatigue  analysis  of  the  thick-section  composite 
structures. 

•  Develop  new  phenomenological  and  micromechanical  fatigue  modeling  capabilities  that 
can  account  for  microstructural  effects  in  the  thick-section  structures  (e.g.,  interphasial 
effects,  interfacial  damage,  stacking  sequence  effects,  etc.). 

6.4  Recommendations  for  Studies  of  Interphasial  Effects  on  LCF.  The  ARL  LCF 

Program  will  help  ARL  engineers  to  meet  the  AAN  goals  if  the  full  potential  of  the  "material-by- 
design"  approach  is  realized.  Section  2.4  includes  some  recommendations  for  research  concerning 
the  effects  of  mechanical  properties  of  the  fiber-matrix  interphases  on  LCF  behavior  of  PMCs 
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and  damage  accumulation  of  composite  structures.  These  recommendations  can  be  summarized 
as  follows. 

•  Develop  experimental  methodology  and  an  appropriate  data  reduction  scheme  for 
characterization  of  the  fiber-matrix  interphase  properties  around  an  individual  fiber  under 
various  strain  rates.  A  new  test  apparatus  called  DILA  and  developed  by  Tanoglu  et  al. 
(to  be  published)  can  be  employed. 

•  Develop  a  physical  understanding  of  the  energy  absorption  mechanisms  activated  in  the 
interphasial  region  between  a  fiber  and  the  matrix  (e.g.,  interfacial  debonding,  interfacial 
void  growth,  fiber  push-out,  and  frictional  fiber  sliding). 

•  Characterize  the  effect  of  energy-absorbing  capabilities  of  the  fiber-matrix  interphase  on 
the  LCF  and  ballistic  performance  of  typical  Army  land  combat  composite  structures. 

•  Develop  models  describing  various  energy-absorbing  mechanisms  activated  in  the  fiber- 
matrix  interphase  and  predict  the  possible  effects  on  the  LCF  fatigue  performance  of 
composites. 

•  Develop  a  physical  understanding  of  the  matrix  toughening  mechanisms  activated  around 
rubber  or  thermoplastic  particles  in  the  matrix  (e.g.,  interfacial  debonding,  interfacial  void 
growth,  plastic  shear  localization,  and  other). 

•  Develop  techniques  for  "interfacial  damage  healing"  (e.g.,  heating  of  magnetic  microscopic 
particles  distributed  around  fibers)  for  different  types  of  PMCs. 

•  Develop  models  describing  various  matrix  toughening  mechanisms  activated  in  around 
particles  in  the  matrix  and  predict  the  possible  effects  on  the  LCF  fatigue  performance  of 
composites. 

6.5  Recommendations  for  Studies  of  Adhesive  Joints.  Section  2.5  includes  some 
recommendations  for  investigation  of  the  role  of  adhesive  joints  in  composites,  which  can  be 
summarized  as  follows. 

•  Identify  the  adhesive  joints  between  dissimilar  composite  materials  which  are  critical  for 
ballistic  and  LCF  performance  of  specific  Army  land  combat  systems.  This  assessment 
should  be  based  on  the  structural  stress  analysis  of  specific  Army  structures. 
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•  Develop  models  describing  the  effects  of  adhesive  joints  on  the  ballistic  damage  tolerance 
of  various  Army  land  combat  systems. 

•  Develop  phenomenological  and  micromechanical  fatigue  models  that  can  describe  the  LCF 
mechanics  of  damaged  hybrid  structures  and  predict  the  LCF  performance  of  damaged 
composite  systems. 

•  Develop  experimental  methodology  and  an  appropriate  data  reduction  scheme  for 
characterization  of  the  energy  absorption  capabilities  of  various  adhesive  joints  under  the 
strain  rates  encountered  during  a  ballistic  impact  or  in  LCF  conditions. 

•  Select  appropriate  test  fixtures  that  can  be  used  for  fatigue  and  fracture  toughness 
characterization  of  the  test  coupons  representing  the  damaged  and  repaired  structural 
elements. 

•  Characterize  different  types  of  fatigue  damage,  damage  evolution  rates,  and  determine  the 
dominant  fatigue  damage  mechanisms  in  the  damaged  and  repaired  composites. 

•  Develop  conceptual  maps  describing  the  accumulation  processes  of  fatigue  damage  in  the 
damaged  and  repaired  composites  under  a  wide  range  of  loading  conditions. 

•  Develop  phenomenological  and  micromechanical  fatigue  models  that  can  be  used  to 
predict  the  stiffness  degradation  rates  and  fatigue  life  of  the  damaged  and  repaired 
composite  systems. 

6.6  Recommendations  for  Investigation  of  Environmental  Effects.  Section  2.6  includes 
some  recommendations  for  future  studies  of  environmental  effects  on  composites,  which  can  be 
summarized  as  follows. 

•  Identify  the  possible  ranges  of  temperatures  and  humidity  levels  and  types  of  chemical 
environment  to  which  typical  Army  land  combat  systems  may  be  exposed  to  during 
service  life. 

•  Identify  typical  groups  of  physical/chemical  conditions  for  various  Army  systems  to 
determine  whether  the  mechanical  fatigue  loading  and  ballistic  impacts  will  be 
complemented  with  thermomechanical  loads,  moisture  effects,  chemical  corrosion, 
mechanical  wear,  etc. 
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•  Develop  experimental  methodology  for  fatigue  testing  of  composite  structures  that  have 
been  exposed  to  various  environmental  conditions. 

•  Develop  experimental  testing  fixtures  for  fatigue  characterization  of  composite  structures 
that  are  being  exposed  to  various  environmental  conditions. 

•  Develop  models  describing  the  effects  of  material  heterogeneity,  re-existing  defects  and 
developing  interfacial  and  matrix  cracking  on  the  absorption  and  propagation  of  moisture. 

•  Develop  models  describing  the  effects  of  material  heterogeneity  such  as  existing  defects 
and  developing  interfacial  and  matrix  cracking  on  the  degradation  of  mechanical  properties 
due  to  the  moist  or  chemically  active  environment. 

•  Perform  experimental  and  theoretical  characterization  of  beneficial  effects  caused  by  the 
use  of  various  protective  coatings  on  the  reduction  of  absorption  and  propagation  of 
moisture. 
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1.  Introduction 


The  MRLife  simulation  code1  is  developed  by  the  Materials  Response  Group  (Virginia 
Polytechnic  Institute)  for  the  simulation  of  performance  and  fatigue  life  prediction  for  composite 
laminates.  The  following  summary  of  the  code's  capabilities  is  based  on  the  MRLife  11  User 
Manual.1  As  a  result,  this  review  is  not  intended  to  reflect  the  current  development  efforts  at 
Virginia  Polytechnic  Institute  and  State  University,  Blacksburg,  VA. 

MRLife  1 1  is  suited  for  fatigue  analysis  of  a  wide  variety  of  problems.  Such  problems  include 
delamination  and  failure  of  notched  and  unnotched  materials  with  or  without  moisture  diffusion. 
The  effects  caused  by  thermal  loads,  creep,  stress  relaxation  and  aging  can  also  be  accounted  for 
in  the  analysis.  The  stacking  sequence  in  a  composite  may  vary  and  include  up  to  28  plies.  This 
capability  allows  one  to  assess  the  fatigue  performance  of  56  plies  in  a  thick  symmetric  laminate. 
The  reduction  in  stiffness  and  strength  of  each  ply  is  governed  separately.  The  degradation  rate 
equations  involve  polynomial  functions  that  may  be  different  for  each  ply 

The  continuum  mechanics  representation  of  stiffness  change  requires  knowledge  of  the 
lamina  phenomenological  constants  that  characterize  the  intralaminar  damage  effects.  The 
continuum  damage  mechanics  parameters  depend  on  the  crack  spacing  evolution  among  other 
things.  A  power  law  approximates  the  crack  density.  Micromechanical  calculations  of  the 
lamina  properties  may  take  into  consideration  the  transversely  isotropic  properties  of  different 
fiber  sizings. 

2.  Micromechanical  Modeling 

2.1  Evaluation  of  Mechanical  Properties.  A  concentric  cylinder  (CC)  model  of  Pagano  and 
Tandon1  is  implemented  to  analyze  the  mechanical  properties  of  a  long  fiber  surrounded  by  a 
sheath  of  matrix  material.  The  boundary  of  the  outer  cylinder  is  subjected  to  average  strains  £°. 

The  composite  stress  field  c %  is  determined  by  the  volume  averaging  of  the  stress  field  over  the 
fiber  and  matrix.  The  effective  elastic  moduli  are  evaluated  by  setting  the  strain  and  expansion 
strain  components  equal  to  zero,  excluding  one  strain  component  each  time.  The  free-edge  effects 
are  not  taken  into  account.  The  Halpin-Tsai  equations  provide  another  approximate  scheme  for 


1  Case,  S.  W.,  and  K.  L.  Reifsnider.  MRLife  11  -  A  Strength  and  Life  Prediction  Code  for  Laminated  Composite 
Materials.  Materials  Response  Group,  Virginia  Polytechnic  Institute  and  State  University,  1998. 
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determining  the  mechanical  properties  of  each  lamina  along  the  fiber  direction  by  using  the  rule- 
of-mixture  approach. 

The  Young's  modulus  in  the  transverse  direction  can  be  evaluated  by  a  number  of  models 
(e.g.,  Chamis's  model,2  Gibson's  model,  the  Halpin-Tsai  equations  and  the  CC  model).  The  shear 
modulus  can  be  estimated  by  the  models  of  Chamis,2  Christensen,3  Gibson,4  Halpin  and  Tsai.5 
The  axial  Poisson's  ratio  is  calculated  by  using  either  the  rule-of-mixtures  or  the  CC  model. 
Thermal  and  environmental  effects  are  taken  into  account  via  the  expansion  strain.  The  axial  and 
transverse  expansion  coefficients,  a  1,2  and  P12  are  evaluated  by  the  CC  model2  and  by  the 
Schapery  model.1  The  effectiveness  of  all  models  is  demonstrated  by  comparing  the  results  with 
the  exact  solution  of  Averill  and  Carman1  for  hexagonally  packed  fibers. 

2.2  The  Tensile  Strength  Models.  The  tensile  strength  of  polymer  matrix  composites 
(PMCs)  in  the  fiber  direction  is  evaluated  by  the  model  developed  by  Gao  and  Reifsnider.1  The 
model  is  based  bn  the  probability  analysis  carried  out  by  Batdorf.1  Batdorf  considered  N  fibers 
surrounded  by  the  matrix  material.  Damage  in  the  composite  system  is  assumed  to  involve  only 
the  fiber  breakage  characterized  by  so-called  singlets,  doublets,  or  i-plets.  The  fiber  failures  are 
assumed  to  conform  to  a  two-parameter  Weibull  representation.  The  probability  of  failure  is 
approximated  by  employing  Reifsnider's  formula  for  the  reliability  of  a  fiber  having  a  linear  stress 
variation.1  As  a  result,  one  may  estimate  the  number  of  i-plets,  Q,  and  construct  a  schematic 
diagram  for  several  i-plets  as  a  function  of  the  applied  stress.  The  envelope  of  intersection 
points  formed  defines  the  set  of  unstable  fiber  breakage  that  lead  to  the  composite  failure.  The 
failure  stress  is  given  by  the  lowest  load  at  which  any  unstable  i-plet  lies  on  the  envelope  formed. 

The  broken  fibers  induce  stress  concentration  and  interfacial  debonding  close  to  the  fracture 
site.  In  the  Gao-Reifsnider  model,1  the  stress  concentrations  and  the  ineffective  lengths  for  each 
group  of  adjacent  fiber  fractures  are  predicted  by  the  shear-lag  theory.  The  core  of  broken 
fiber(s)  is  flanked  by  a  layer  of  unbroken  fibers  and  the  outer  layer  of  a  homogenized  material 
with  effective  properties.  It  is  assumed  that  the  core  of  i  broken  fibers  is  approximated  by  a 
circular  cross  section  with  the  Young's  modulus  determined  by  the  rule-of-mixtures.  The  degree 


1  Case,  S.  W.,  and  K.  L.  Reifsnider.  MRLifell  -  A  Strength  and  Life  Prediction  Code  for  Laminated  Composite 
Materials.  Materials  Response  Group,  Virginia  Polytechnic  Institute  and  State  University,  1999. 

2  Chamis,  C.  C.  "Simplified  Composite  Micromechanics  Equations  for  Strength,  Fracture  Toughness,  Impact 
Resistance  and  Environmental  Effects."  NASA  Technical  Memorandum  83696,  NASA  Lewis  Research  Center, 
Cleveland,  OH,  1984. 

3  Christensen,  R.  M.  Effective  Modulii  of  Cylindrical  and  Lamellar  Systems.  Mechanics  of  Composite  Materials, 
pp.  73 1 05.  New  York,  NY:  Wiley,  1 979. 

4  Gibson,  R.  F.  Principles  of  Composite  Material  Mechanics.  New  York,  NY:  McGraw  Hill,  1994. 

5  Jones,  R.  M.  Mechanics  of  Composite  Materials.  New  York,  NY:  Hemisphere  Publishing  Co.,  1975. 
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of  interfacial  debonding  is  defined  by  a  shear  parameter,  which  can  represent  complete  fiber- 
matrix  debonding  and  elastic-perfectly-plastic  behavior.  The  ineffective  length  is  defined  as  the 
distance  required  for  the  stress  on  the  adjacent  fibers  to  reach  95%  of  the  far-field  stress. 

The  tensile  strength  of  PMCs  in  the  fiber  direction  can  be  also  evaluated  by  Reifsnider's 
model.1  This  model  is  also  based  on  the  Batdorfs  analysis  similarly  to  the  Gao-Reifsnider 
model;1  however,  the  failure  criterion  is  modified.  Another  modification  of  the  Gao-Reifsnider 
model  involves  some  changes  in  the  approximate  geometry  calculations  and  the  effective  stiffness 
of  the  core  of  broken  fibers.1  The  tensile  strength  of  ceramic  matrix  composites  in  the  fiber 
direction  is  evaluated  by  the  model  developed  by  Curtin.1  In  this  model,  a  ceramic  matrix  is 
reinforced  with  uniaxial  fibers  conforming  to  the  two-parameter  Weibull  function.  The  initial 
elastic  properties  of  the  composite  are  estimated  by  the  rule  of  mixtures.  The  load-bearing 
capacity  of  the  composite  is  adjusted  after  matrix  cracking  and  individual  fiber  breakage, 
depending  on  the  fiber  slip  length. 

2.3  The  Compressive  Strength  Models.  The  compressive  strength  of  PMCs  in  the  fiber 
direction  is  evaluated  by  the  model  developed  by  Xu  and  Reifsnider.1  It  employs  a  beam-on- 
elastic  foundation  model  to  consider  the  critical  load  for  fiber  microbuckling.  Interfacial  slipping 
of  the  matrix  is  taken  into  account  in  this  model.  The  half-wavelengths  of  buckled  fibers  are 
determined  by  applying  the  minimum  fiber  buckling  load  condition.  The  cylindrical  fibers  are 
approximated  by  the  square  beams  having  the  same  cross-sectional  area.  The  stiffness  of  the 
foundation  is  determined  through  an  elasticity  solution  to  a  foundation  model  problem. 

The  compressive  strength  of  PMCs  in  the  fiber  direction  can  be  also  determined  by  the  model 
developed  by  Fleck  and  Budiansky.1  This  model  takes  into  account  the  plastic  deformation  via 
kink  bands.  The  Ramburg-Osgood  relation  is  employed  to  model  the  material  nonlinearity.  In 
this  model,  effects  of  the  initial  misalignment  and  variations  in  the  shear  angle  are  accounted  for. 
The  limiting  load  for  fiber  crushing  can  be  estimated  by  the  rule  of  mixtures. 

2.4  The  Transverse  and  Shear  Strength  Models.  The  transverse  tensile  strength  can  be 
estimated  by  the  model  introduced  by  Gibson.2  The  model  is  based  on  the  standard  rule  of 
mixtures,  so  it  does  not  take  into  account  the  changes  in  fiber  packing.  More  complex  models,2 
which  are  also  implemented,  can  simulate  the  effects  of  hexagonal  and  square  packing.  The  in¬ 
plane  shear  strength  can  be  evaluated  by  the  three  analogous  models  suggested  by  Gibson.2  The 


1  Case,  S.  W.,  and  K.  L.  Reifsnider.  MRLifell  -  A  Strength  and  Life  Prediction  Code  for  Laminated  Composite 
Materials.  Materials  Response  Group,  Virginia  Polytechnic  Institute  and  State  University,  1998 

2  Gibson,  R.  F.  Principles  of  Composite  Material  Mechanics.  New  York,  NY:  McGraw  Hill,  1994. 
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Gibson's  family  of  models  also  provides  an  estimate  for  the  transverse  compression  strength.  It 
is  also  based  on  the  rule  of  mixtures  approach. 

3.  Methodology  of  Damage  Analysis 

In  the  damage  models  implemented,  several  important  features  of  the  complex  physical 
behavior  of  composite  systems  are  taken  into  account.  The  damage  modeling  is  carried  out  for 
various  fatigue  loadings,  including  the  influence  of  creep  and  aging.  The  combined  effects  of  the 
loading  conditions  considered  are  characterized  by  using  a  set  of  damage  accumulation  concepts. 

The  transverse  matrix  cracking  in  laminate  composites  is  described  by  Talreja's  tensorial 
representation1  described  in  section  5.1.  In  MRLifel  1,  delamination  is  assumed  to  be  caused  by 
the  interlaminar  stresses  at  a  free  edge.  The  strain  energy  release  rate,  G,  for  the  interlaminar 
crack  is  estimated  by  a  fracture  mechanics  approach  in  conjunction  with  the  laminated  plate 
theory.  The  estimate  is  independent  of  the  delamination  length,  but  it  takes  into  account  the 
thickness  of  the  laminate.  The  O'Brien-Paris  law1  is  employed  for  the  interlaminar  crack  growth 
from  a  delaminated  edge  under  fatigue  loading.  This  model  for  the  interlaminar  fatigue  crack 
growth  is  similar  to  the  Paris  law  for  fatigue  crack  growth  in  metals. 

To  model  accumulation  of  damage,  MRLifel  1  employs  a  scheme  developed  by  Reifsnider  et 
al.1  for  composites  under  fatigue  loading.  It  is  first  postulated  that  the  damage  can  be 
characterized  by  the  remaining  strength,  which  is  a  function  of  the  level  of  load  and  a  generalized 
time.  The  equivalence  between  different  fractions  of  fatigue  life,  which  correspond  to  the  same 
reduction  in  remaining  strength  under  different  loads,  is  also  postulated.  The  remaining  fatigue 
life  at  the  load  applied  is  determined  by  the  amount  of  generalized  time  required  to  reduce  the 
remaining  strength  to  the  applied  load  level.  As  a  result,  the  effect  of  changes  in  loading  may  be 
taken  into  account  by  adding  the  respective  reductions  in  remaining  strength.  Since  the  strength 
reduction  curves  may  be  nonlinear,  the  remaining  strength  and  corresponding  life  prediction 
calculations  are  path  dependent. 

MRLifel  1  has  an  extensive  library  of  failure  criteria.  Once  an  appropriate  criterion  for  failure 
is  chosen,  the  normalized  remaining  strength  can  be  defined  as  an  internal  state  variable  for  a 
damaged  material  system.  The  Kachanov's  continuity  function  \\t  is  the  second  state  variable 
related  to  the  Helmholtz  free  energy.  A  specific  damage  accumulation  process  for  a  particular 


'  Case,  S.  W.,  and  K.  L.  Reifsnider.  MRLifel  1  -  A  Strength  and  Life  Prediction  Code  for  Laminated  Composite 
Materials.  Materials  Response  Group,  Virginia  Polytechnic  Institute  and  State  University,  1998 
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failure  mode  has  its  own  damage  kinetics,  which  is  described  by  special  rate  equations  (e.g.,  a 
power  law,  etc). 

The  damage  kinetic  equations  require  inputs  about  the  evolution  of  damage  and  its  effect  on 
the  failure  of  representative  material  elements.  The  concepts  of  representative  “critical”  and 
“subcritical”  material  elements  are  used  to  characterize  the  accumulation  of  damage  within  the 
composite  laminate.  The  0°  plies  in  a  cross-ply  laminate  represent  an  example  of  “critical” 
elements,  while  90°  plies  represent  “subcritical”  elements  on  the  ply  level.  The  micromechanical 
effects  can  be  also  incorporated  into  such  framework. 
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M  LUCIANO 
ELOGSDEN 
T  LOUZEIRO 

PICATINNY  ARSENAL  NJ 
07806-5000 

4  COMMANDER 

US  ARMY  ARDEC 
AMSTAARCC 
G PAYNE 
J  GEHBAUER 
CBAUUEU 
HOPAT 

PICATINNY  ARSENAL  NJ 
07806-5000 


66 


NO.  OF 

COPIES  ORGANIZATION 


1  COMMANDER 

US  ARMY  ARDEC 
AMSTAARCCHP 
JLUTZ 

PICAT1NNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 

AMSTAARFSFT 

C  UVECCHIA 

PIC  ATINNY  ARSENAL  NJ 

07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARQACTC 
C PATEL 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  M 

D  DEMELLA 
FDIORIO 

PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARFSA 
A  WARN  ASH 
BMACHAK 
MCHIEFA 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARFSPG 
M  SCHIKSNIS 
D  CARLUCCI 
PICATINNY  ARSENAL  NJ 
07806-5000 


NO.  OF 

COPIES  ORGANIZATION 


1  COMMANDER 
US  ARMY  ARDEC 
AMSTAARFSPA 
PKISATSKY 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARCCHC 
H  CHAN  IN 

S  CHICO 

PICATINNY  ARSENAL  NJ 
07806-5000 

9  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARCCHB 
P  DONADIA 
FDONLON 
P  VALENTI 
C KNUTSON 
GEUSTICE 
S PATEL 
G  WAGNECZ 
RSAYER 
F CHANG 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  CCL 
FPUZYCKI 
R  MCHUGH 
D  CONWAY 
EJAROSZEWSKI 
R  SCHLENNER 
MCLUNE 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARQACT 
DRIGOGUOSO 
PICATINNY  ARSENAL  NJ 
07806-5000 
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1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  SRE 
DYEE 

PIC  ATINN  Y  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  WET 
TSACHAR 
BLDG  172 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
SMCAR  ASF 

PICATINNY  ARSENAL  NJ 
07806-5000 

I  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARWELF 
INTELLIGENCE  SPECIALIST 
M  GUERRIERE 
PICATINNY  ARSENAL  NJ 
07806-5000 

II  PROJECT  MANAGER 
TANK  MAIN 
ARMAMENT  SYSTEMS 
SFAEGSSCTMA 

R  MORRIS 
CKIMKER 
D  GUZOWICZ 
EKOPACZ 
RROESER 
R DARCY 
R  MCDANOLDS 
LDULISSE 
C  ROLLER 
J  MCGREEN 
B  PATTER 

PICATINNY  ARSENAL  NJ 
07806-5000 


NO.  OF 

COPIES  ORGANIZATION 


2  PEO  FIELD  ARTILLERY  SYSTEMS 
SFAEFASPM 

H  GOLDMAN 
T  MCWILLIAMS 
PICATINNY  ARSENAL  NJ 
07806-5000 

6  PM  SAD ARM 

SFAEGCSS  SD 
COL  B  ELLIS 
M  DEVINE 
R  KOWALSKI 
W  DEMASSI 
J  PRITCHARD 
SHROWNAK 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
PRODUCTION  BASE 
MODERN  ACTY 
AMSMC  PBM  K 
PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 
USARMYTACOM 

PM  TACTICAL  VEHICLES 
SFAETVL 
SFAETVM 
SFAETVH 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

USARMYTACOM 
PM  ABRAMS 
SFAE  ASM  AB 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

USARMYTACOM 

PMBFVS 

SFAE  ASM  BV 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 
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COMMANDER 

1 

COMMANDER 

US  ARMYTACOM 

US  ARMYTACOM 

PMAFAS 

PM  GROUND  SYSTEMS 

SFAE  ASM  AF 

INTEGRATION 

6501  ELEVEN  MILE  RD 

SFAEGCSS  W GSI 

WARREN  MI  48397-5000 

RLABATILLE 

6501  ELEVEN  MILE  RD 

COMMANDER 

U  S  ARMY  TACOM 

WARREN  MI  48397-5000 

PM  SURV  SYS 

1 

COMMANDER 

SFAE  ASM  SS 

US  ARMY  TACOM 

TDEAN 

CHIEF  ABRAMS  TESTING 

SFAEGCSS  WGSIM 

SFAEGCSS  WAB  QT 

D COCHRAN 

T  KRASK3EWICZ 

6501  ELEVEN  MILE  RD 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

WARREN  MI  48397-5000 

COMMANDER 

1 

COMMANDER 

US  ARMYTACOM 

US  ARMYTACOM 

PMRDT&E 

AMSTA  SF 

SFAEGCSS  WAB 

JGODELL 

WARREN  MI  48397-5000 

6501  ELEVEN  MILE  RD 

1 

COMMANDER 

WARREN  MI  48397-5000 

SMCWV  QAEQ 

B  VANINA 

COMMANDER 

BLDG  44 

US  ARMYTACOM 

WATER  VLIET  ARSENAL 

PM  SURVIVABLE  SYSTEMS 

SFAE  GCSS  W  GSI H 

WATERVLIETNY  12189-4050 

M  RYZYI 

14 

COMMANDER 

6501  ELEVEN  MILE  RD 

US  ARMYTACOM 

WARREN  MI  48397-5000 

ASMTATRR 

J  CHAPIN 

COMMANDER 

R  MCCLELLAND 

US  ARMYTACOM 

D  THOMAS 

PMBFV 

J BENNETT 

SFAE  GCSS  WBV 

D  HANSEN 

S  DAVIS 

AMSTA  JSK 

6501  ELEVEN  MILE  RD 

S GOODMAN 

WARREN  MI  48397-5000 

J  FLORENCE 

KIYER 

COMMANDER 

J  THOMSON 

U  S  ARMY  TACOM 

AMSTA  TRD 

PM  LIGHT  TACTICAL 

D  OSTBERG 

VEHICLES 

L  HINOJOSA 

AMSTATRS 

BRAJU 

AJJ  MILLS  MS  209 

AMSTA  CSSF 

6501  ELEVEN  MILE  RD 

H  HUTCHINSON 

WARREN  MI  48397-5000 

F  SCHWARZ 

WARREN  MI  48397-5000 
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1  COMMANDER 
SMCWV  SPM 
T  MCCLOSKEY 
BLDG  253 

WATER  VLIET  ARSENAL 
WATER  VUET  NY  12189-4050 

1 0  BENET  LABORATORIES 

AMSTAARCCB 
RFISCELLA 
GD  ANDREA 
M  SCAVULO 
G  SPENCER 
P  WHEELER 
K  MINER 
J  VASILAKIS 
G  FRIAR 
RHASENBEIN 
SMCAR  CCB  R 
S  SOPOK 

WATERVUETNY  12189 

2  TSM  ABRAMS 
ATZKTS 
SJABURG 

W  MEINSHAUSEN 
FT  KNOX  KY  40121 

3  ARMOR  SCHOOL 
ATTN  ATZKTD 
RBAUEN 
JBERG 
APOMEY 

FT  KNOX  KY  40121 

2  HQ  IOC  TANK  AMMO  TEAM 

AMSIO  SMT 
R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 

1  DIRECTOR 

U  S  ARMY  AMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BUILDING  5300 

REDSTONE  ARSENAL  AL  35898 


4  DIRECTOR 

USARMYCECOM 
NIGHT  VISION  AND  ELECTRONIC 
SENSORS  DIRECTORATE 
AMSEL  RD  NV  CM  CCD 
R  ADAMS 
R  MCLEAN 
A  YINGST 
AMSEL  RDNVVISP 
E  JACOBS 
10221  BURBECK  RD 
FT  BELVOIR  VA  22060-5806 

2  CDR  USA  AMCOM 

AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604-5577 

1  US  ARMYCRREL 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  US  ARMYCERL 
RLAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  61 822 

2  US  ARMY  CORP  OF  ENGINEERS 
CERD  C  T  LIU 
CEWETTTAN 

20  MASS  AVENW 
WASHINGTON  DC  20314 

10  DIRECTOR 

US  ARMY  NATL  GRND  INTEL  CTR 

DLEITER 

SEITELMAN 

MHOLTUS 

M  WOLFE 

S  MINGLEDORF 

H  C  ARDLEIGH 

J  GASTON 

W  GSTATTENBAUER 
R  WARNER 
J  CRIDER 

220  SEVENTH  STREET  NE 
CHARLOTTESVILLE  VA  22091 
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6  US  ARMY  SBCCOM 

SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
J  WARD 

MARINE  CORPS  TEAM 

J  MACKIEWICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

SSCNC  WST 

WNYKVIST 

T  MERRILL 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

1  US  ARMY  COLD 

REGIONS  RSCH  & 

ENGRNG  LAB 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  SYSTEM  MANAGER  ABRAMS 

BLDG  1002  RM 110 
ATZK  TS  LTC  J  H  NUNN 
FT  KNOX  KY  40121 

9  US  ARMY  RESEARCH  OFFICE 

ACROWSON 
J CHANDRA 
H EVERETT 
J PRATER 
R  SINGLETON 
G ANDERSON 
D  STEPP 
D  KISEROW 
J CHANG 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  DIRECTORATE  OF  CMBT 

DEVELOPMENT 
CKJORO 

320  ENGINEER  LOOP  STE  141 
FT  LEONARD  WOOD  MO  65473-8929 


1  COMMANDANT 

U  S  ARMY  FIELD  ARTILLERY  CTR 
AT  FT  SILL 

ATFS  CD  LTC  BUMGARNER 
FT  SILL  OK  73503  5600 

1  CHIEF  US  AIC 

LTC  TJ  CUMMINGS 
ATZBCOM 

FT  BENNING  GA  31905-5800 

1  NAVAL  AIR  SYSTEMS  CMD 
J  THOMPSON 
48142  SHAW  RD  UNIT  5 
PATUXENT  RIVER  MD  20670 

1  NAVAL  SURFACE  WARFARE  CTR 

DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 

1  NAVAL  SURFACE  WARFARE  CTR 

TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

3  NAVAL  RESEARCH  LAB 
I WOLOCK  CODE  6383 
R  BADAIIANCE  CODE  6304 
LGAUSE 

WASHINGTON  DC  20375 

1  NAVAL  SURFACE  WARFARE  CTR 
CRANE  DIVISION 

M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CADEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFTELD  CODE  1730 
BETHESDA  MD  20084 

2  NAVAL  SURFACE  WARFARE  CTR 

U  SORATHIA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 
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DAVID  TAYLOR  RESEARCH  CTR 
SHIP  STRUCTURES  &  PROTECTION 
DEPARTMENT  CODE  1702 
J  CORRADO 
BETHESDA  MD  20084 

DAVID  TAYLOR  RESEARCH  CTR 
R  ROCKWELL 
WPHYELLAIER 
BETHESDA  MD  20054-5000 

OFFICE  OF  NAVAL  RESEARCH 
D  SIEGEL  CODE  351 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 

NAVAL  SURFACE  WARFARE  CTR 
J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
RD  COOPER  CODE  G32 
JFRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 

NAVAL  SEA  SYSTEMS  CMD 
DUESE 

2531  JEFFERSON  DAVIS  HIGHWAY 
ARLINGTON  VA  22242-5160 

NAVAL  SURFACE  WARFARE 
M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

OFFICE  OF  NAVAL  RES 
J  KELLY 

800  NORTH  QUINCEY  ST 
ARLINGTON  VA  22217-5000 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3  A  LEGGETT  CIR 
BETHESDA  MD  20054-5000 


1  NAVSEA  OJRI 

PEO  DD21  PMS500 
G  CAMPONESCHI 
2351  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5165 

1  EXPEDITIONARY  WARFARE 

DIVN85 
FSHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  AFRLMLBC 

2941  P  STREET  RM 136 
WRIGHT  PATTERSON  AFB  OH 
45433-7750 

1  AFRLMLSS 
R  THOMSON 

2179  12TH  STREET  RM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F ABRAMS 
J  BROWN 
BLDG  653 

2977  P  STREET  STE  6 
WRIGHT  PATTERSON  AFB  OH 
45433-7739 

1  AFRL  MLS  OL 

7278  4TH  STREET 
BLDG  100  BAY  D 
L  COULTER 

HILL  AFB  UT  84056-5205 
1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD  R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

1  DEFENSE  NUCLEAR  AGENCY 

INNOVATIVE  CONCEPTS  DIV 
RROHR 

6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 
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WATERWAYS  EXPERIMENT 
D  SCOTT 

3909  HALLS  FERRY  RDSCC 
VICKSBURG  MS  39180 

DARPA 

M  VANFOSSEN 
SWAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

SERDP  PROGRAM  OFC  PM  P2 
C  PELLERIN 
B  SMITH 

901  N  STUART  ST  SUITE  303 
ARLINGTON  VA  22203 

FAA 

MIL  HDBK 17  CHAIR 
LILCEWICZ 
1601  UND  AVE  SW 
ANM115N 
RENTON  VA  98055 

FAA  TECH  CENTER 
D  OPUNGER  AAR  431 
P  SHYPRYKEVICH  AAR  431 
ATLANTIC  CITY  NJ  08405 

OFC  OF  ENVIRONMENTAL  MGMT 
US  DEPT  OF  ENERGY 
P  RJTZCOVAN 
19901  GERMANTOWN  RD 
GERMANTOWN  MD  20874-1928 

LOS  ALAMOS  NATL  LAB 
F  ADDESSIO 
MS  B216 
PO  BOX  1633 
LOS  ALAMOS  NM  87545 

OAK  RIDGE  NATL  LAB 
RM  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 


5  DIRECTOR 

LAWRENCE  LIVERMORE 
NATL LAB 
R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

7  NIST 

RPARNAS 
J  DUNKERS 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
D  HUNSTON  MS  8543 
J  MARTIN  MS  8621 
D  DUTHINH  MS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

1  OAK  RIDGE  NATL  LAB 

C  EBERLE  MS  8048 
PO  BOX  2009 
OAK  RIDGE  TN  37831 

1  OAK  RIDGE  NATL  LAB 

CD  WARREN  MS  8039 
PO  BOX  2009 
OAK  RIDGE  TN  37922 

4  DIRECTOR 

SANDIA  NATL  LABS 
APPLIED  MECHANICS  DEPT 
DIVISION  8241 
WKAWAHARA 
KPERANO 
D  DAWSON 
PN1ELAN 
PO  BOX  969 

LIVERMORE  CA  94550-0096 

1  LAWRENCE  LIVERMORE 
NATIONAL  LAB 
M  MURPHY 
POBOX  808  L  282 
LIVERMORE  CA  94550 
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NASA  LANGLEY  RESEARCH  CTR 

MS  266 

AMSRLVS 

WELBER 

F  BARTLETT  JR 

G FARLEY 

HAMPTON  VA  23681-0001 

NASA  LANGLEY  RESEARCH  CTR 
T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

USDOT  FEDERAL  RAILROAD 
RDV31  M  FATEH 
WASHINGTON  DC  20590 

DOTFHWA 

JSCALZI 

400  SEVENTH  ST  SW 
3203  HNG  32 
WASHINGTON  DC  20590 

FHWA 

EMUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 

CENTRAL  INTELLIGENCE  AGENCY 

OTIWDAGGT 

WL  WALTMAN 

PO  BOX  1925 

WASHINGTON  DC  20505 

MARINE  CORPS  INTEL  ACTY 
DKOSITZKE 

3300  RUSSELL  RD  SUITE  250 
QUANTICO  VA  22134-501 1 

NATL  GRND INTELUGENCE  CTR 

DIRECTOR 

IANGTMT 

220  SEVENTH  ST  NE 

CHARLOTTESVILLE  VA  22902-5396 

DIRECTOR 

DEFENSE  INTELUGENCE  AGENCY 
TA  5  K  CRELUNG 
WASHINGTON  DC  20310 


1  GRAPHITE  MASTERS  INC 

J  WILLIS 

3815  MEDFORD  ST 

LOS  ANGELES  CA  90063-1900 

1  ADVANCED  GLASS  FIBER  YARNS 

T  COLLINS 

281  SPRING  RUN  LNSTEA 
DOWNINGTON  PA  19335 

1  COMPOSITE  MATERIALS  INC 

D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  COMPOSITE  MATERIALS  INC 

R  HOLLAND 
11  JEWEL  COURT 
OR1NDACA  94563 

1  COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

2  COMPOSIX 
D BLAKE 

L  DIXON 

120  0  NEILL  DR 

HEBRUN  OHIO  43025 

4  CYTEC  FIBERITE 

R DUNNE 
DKOHU 
MGILUO 
RMAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

2  SIMULA 

J  COLTMAN 
RHUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 

1  SIOUX  MFG 

BKRIEL 
PO  BOX  400 
FT  TOTTEN  ND  58335 
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2  PROTECTION  MATERIALS  INC 
M  MILLER 

F  CRILI-EY 
14000  NW  58  CT 
MIAMI  LAKES  EL  33014 

3  FOSTER  MILLER 
J  J  GASSNER 

M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1 196 

1  ROM  DEVELOPMENT  CORP 
ROMEARA 

136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI 02840 

2  TEXTRON  SYSTEMS 
T FOLTZ 

M TREASURE 
201  LOWELL  ST 
WILMINGTON  MA  08870-2941 

1  JPS  GLASS 

L CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 

9113  LESAINT  DR 
FAIRFIELD  OH  45014 

2  MHUKEN  RESEARCH  CORP 
H  KUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 

J  SANTOS 
PO  BOX  1425 
COVENTRY  RI  02816 


1  BATTELLE 

CR  HARGREAVES 
505  KING  AVE 
COLUMBUS  OH  43201-2681 

2  BATTELLE  NATICK  OPERATIONS 
J CONNORS 

BHALPIN 

209  W  CENTRAL  ST 
STE  302 

NATICK  MA  01760 

1  BATTELLE  NW  DOE  PNNL 

T  HALL  MS  K231 
BATTELLE  BLVD 
RICHLAND  WA  99352 

3  P  ACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
RICHLAND  WA  99352 

1  ARMTEC  DEFENSE  PRODUCTS 
S  DYER 

85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

2  ADVANCED  COMPOSITE 
MATLS  CORP 
PHOOD 

J  RHODES 

1525  S  BUNCOMBE  RD 
GREER  SC  29651-9208 

2  GLCC  INC 
J  RAY 

M  BRADLEY 

103  TRADE  ZONE  DR 

STE26C 

WEST  COLUMBIA  SC  29170 

2  AMOCO  PERFORMANCE 
PRODUCTS 
M  MICHNO  JR 
J  BANISAUKAS 
4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 


75 


NO.  OF 

COPIES  ORGANIZATION 


1  SAIC 

M  PALMER 

2109  AIR  PARK  RD  S  E 

ALBUQUERQUE  NM  87106 

1  SAIC 

ATTN  G  CHRYSSOMALLIS 
3800  W  80TH  ST  STE  1090 
BLOOMINGTON  MN  55431 

1  AAI  CORPORATION 

DR  T  G  STASTNY 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  JOHN  HEBERT 

PO  BOX  1072 

HUNT  VALLEY  MD  21030-0126 

1 2  ALLIANT  TECHS YSTEMS  INC 

C  CANDLAND 
CAAKHUS 
R BECKER 
B  SEE 

N  VLAHAKUS 
RDOHRN 
S  HAGLUND 
D  FISHER 
W  WORRELL 
R  COPENHAFER 
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5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to  organization, 
technical  content,  format,  etc.) _ _ _ _ _ 


Organization 

CURRENT  Name  E-mail  Name 

ADDRESS  _ _ _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


